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The  goal  of  this  research  was  to  increase  available  signal-to- 
noise  ratio  (SNR)  in  magnetic  resonance  imaging  (MRI)  by  applying 
specific   knowledge   of  the   imaging   system   to   improve   receiver  probes 
(coils)   and  receiving  hardware.      A  brief  history  of  improvements   in 
MRI  receiver  and  coil   design  is  presented,   including  the  transition 
from  large     linear  volume  coils  to  local  surface  coils  and  quadrature 
volume  coils.      Then   quadrature   surface  coils   are   introduced   and 
finally   multi-coil   arrays  with  independent  acquisition   systems.      The 
research   covers   improvements   in   these   areas   and   begins   with   a 
surface  coil  which  is  adjustable  in  size  to  optimize  performance  given 
the  region  of  interest.     By  careful  design  of  trombone-like  coil 
elements,    physical    adjustment   can   be   made   without   electrical 
adjustment.      Second,   new   understanding   of  noise   correlation   and 
crosstalk  between   coils   is   developed   and   applied   to   muli-coil   arrays. 
This  provides  the  ability  to  increase  available   SNR  for  such  systems. 
Third,   a  method  for  optimally  combining  multiple  coils  in  a 
transverse    (extending    perpendicular    to    the    static    magnetic    field) 
array  into  a  single  channel  by  proper  signal  combination  is 


presented.      This  method  is   termed  generalized  quadrature   because  of 
the   similarity   of  the   method   to    standard   quadrature   combination,    but 
with  freedom  in  weighting  and  phasing  in  the  combination  process. 
Fourth,   several  methods  of  manipulating  the  multiple  signals  from   an 
array   to   allow   separation   after   acquisition   are   presented.      These 
methods    require    new    hardware    demands    but    allow    significant 
improvments   in   SNR  for  either  transverse  or  longitudinal   arrays. 
Fifth,   several  novel  design  methods   are  demonstrated,   including   an 
algorithm    for    impedance    matching,    a   generalized    quadrature 
combination   method,   transmission   synchronized   rf  shielding   and   a 
bird-cage   surface   coil.      Finally,   the   potential   future   applications   and 
benefits    of  this   research    are   presented. 
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CHAPTER  1 
INTRODUCTION 

The  Phenomenon  of  Nuclear  Magnetic  Resonance  (NMR) 

The  properties  of  the  nuclei  of  atoms  are  described  by 

quantum  mechanics.     The  significant  property  for  this  work  is  the 

property  of  spin.     Some  nuclei  possess  an  intrinsic  angular 

momentum  or  spin.   Denoting   spin   S,  we  can  describe  its  quantized 

values  as     S  =    — Js(s  +  l)   where  h  is  Planck's  constant  and  s  is  the 
271 

spin  quantum  number  and  is  integral  or  half  integral.      A   second 

quantum   number   is   introduced   to   indicate   as   much   as   allowable 

about  the  direction  of  the  spin  in  a  magnetic  field.     The  convention  is 

h 
to  define  a  magnetic  field  in  the  z-direction  with   S,  =—  mv   being  the 

271 

projection  of  the  spin  onto  the  z-axis.       The  quantum  number  ms  may 

take  on  the  values  of  (s,s-l,s-2,  ...,-s).     A  magnetic  moment  can  be 

thought  of  as  being   the  result  of  an  electric  charge  distribution 

within   a   spinning   nucleus. 

Since  the   spin  magnetic  moment  is  quantized  in   a  magnetic 

field  it  can   be  identified  with  certain  energy   states.     Hydrogen  nuclei 

are   almost  always   the   nuclei   of  interest  in   magnetic   resonance 

imaging   (MRI)   and   their  spin  properties   are  characterized   as   spin    1/2. 


These  nuclei  have  2  allowable  quantum  states  in  a  magnetic  field. 

One   is   higher  energy,   with   the   magnetic   moment   anti-parallel   with 

the  field  and  the  other  of  lower  energy  which  is  parallel.     The  energy 

difference  is   A£  =  — Y^o  where  Y  is  the  gyromagnetic  ratio  and  is 
2n 

defined   empirically   from   this   expression   and   thus   relates   the   angular 
frequency  of  radiation  for  a  transition  between  states,   to  B0,    the 
static   magnetic  induction.     The  population  difference  of  the  two 
energy   states  defines   the   amplitude   of  the  bulk  magnetic  moment 
which   is   available   for  perturbation   and   is   described   by   Boltzman's 

AT 

equation   as    n„-np  =N where  N  is  the  number  of  nuclei,  k  is 

Boltzman's  constant  and  T  is  the  absolute  temperature  of  the  sample. 
The  population  difference,    -    "" ,  is  typically  on  the  order  of  10" 5 

for  typical  MRI  situations  which  is  the  root  of  the  problem  of 
sensitivity  in   MRI. 

In  quantum  mechanics,   a   Hamiltonian   which   describes   the 
energy   states  of  the  system  is   created  in  analogy  with  classical 
mechanics.     The  energy  states  of  a  nuclei  are  in  general  extremely 
complex.     In  the  following  description  all  terms  except  that  of  the 
energy  of  a  magnetic  moment  in  a  magnetic  field  will  be  ignored. 
This  term  is 

H  =  -MH, 
where   H   is  the  spin  Hamiltonian,  M  is  the  magnetic  moment  of  the 
spin  (proportional  to  the  spin  operator  S)  and  H  is  the  magnetic  field. 
By  neglecting   orbital   angular  momentum,   the   behavior  of  the  spin   in 
a  magnetic  field  can  be  readily  observed.     In  the  Heisenberg 
represention   the   equation   of  motion   is 


,Asi(r)  =  [Si(r),H]. 


2% 

This  can  be  rewritten  using  the  commutation  relations  for  the 
components  of  the  spin  operator  and  the  definition  of  the  spin 
Hamiltonian   above   as 

S(O  =  co0S(f)xH(O. 
This  equation  is  an  operator  description  of  the  classical  expression 
which   equates   the   time   derivative   of  angular   momentum    and   the 
applied  torque.     In  NMR  this  expression  is  called  the  Bloch  equation 
and  describes  the  macroscopic  motion  of  the  bulk  magnetic  moment 
in  a  material.     If  the  magnetic  field  is  constant  in  time,  it  is  easily 
shown  that  the  spin  moment  will  precess  or  rotate  about  the  axis  of 
the  applied  field.     The  rate  at  which  this  occurs  is  called  the  Larmor 
frequency.      As   with   classical   systems   which  have   a  natural 
frequency,   excitation   of  the   system   is   most  easily   produced  by 
coupling  to  the  system  at  that  frequency.     In  this  case  a  magnetic 
field  which   varies   at  the   Larmor  frequency   will   strongly  couple   to 
the  system.     It  should  also  be  clear  that  for  a  given  amount  of  power 
used   in   producing   a  time-varying   magnetic   field   at  the  frequency   of 
interest,   a  rotating   magnetic   field  will   most  strongly   interact  with   the 
spin   magnetic   moment.      Therefore,   if  other  parameters   are   held 
constant   a   magnetic   field   described   by   an   exponential,   with   exponent 
i'COor,   will   couple  most  efficiently   and  optimally   produce   the   dynamics 
of  interest.     This  will  be  further  expounded  upon  in  the  following 
sections.      Figure    1-1    demonstrates   the   relationship   of  the   magnetic 
moment   precession,    the   perturbing    rf  magnetic   induction   B]    and  the 
main  fixed  magnetic  induction  B0. 


NMR 


Figure   1-1.     Pictorial  representation  of  the  precession  of  the  spin 
magnetic  moment  and  the  spatial  relationship  with  the  static  field 
Bo  and  the  RF  field  Bl 

In  reality,   as  with   any  oscillator  (other  than  a  perfect 
superconductor),  the  oscillations  will  decrease  with  time.     This  has 
been  termed  relaxation  and  occurs  due  to  a  number  of  effects.     The 
energy   in   the  excited   state   nuclei   which   form  the   net  magnetization 
will  be  dissipated  in  some  way.     One  method  is  by  local  interactions 
of  the  spin  with  its  surroundings  or  lattice.     This  form  of  dissipation 
is   termed   spin-lattice   relaxation   and   the   time   constant  of  this   process 
is  called  Ti.     In  general  terms  each  spin  is  effected  by  other  magnetic 
nuclei  which  fluctuate  to  some  degree  at  the  Larmor  frequency.      At 
typical  MR  field  strengths,   the  rotational  rates  of  medium   sized 
molecules    are   approximately   the   inverse   of   the   Larmor   frequency 


thus  protons  in  this  environment  have  the   most  efficient  relaxation 
and  therefore  the  shortest  Ti   (1). 

The  other  common  relaxation  method,   termed  T2   relaxation 
describes   the  case   where  energy  is  transferred  from  an  excited 
nucleus  to  a  ground  state  nucleus.     If  all  nuclei  were  precessing  at 
the   same   rate,   the   magnetization   would   not  decay   (ignoring   spin- 
lattice  relaxation).     However,  in  real  magnets  and  with  real  samples 
there  is  always  some  inhomogeneity  in  the  magnetic  field,  which  in 
turn   implies   a  difference   in   the   precessional   frequencies   of  various 
spins.     The  net  effect  is  a  loss  of  coherence  which  eventually  leads  to 
no  net  transverse  magnetization.     Again  the  local  magnetic  fields  of 
other  nuclei  are  most  significant  in  this  process.     However,   here   the 
perturbation  of  frequency  is  most  pronounced  for  large,   slow  moving 
molecules   which   therefore   produce   short  T2.     For  typical  samples 
used  in  MRI,  too  little  signal  is  absorbed  by  the  receiver  to  damp 
oscillations,    although   this   phenomenon   has   been   observed.    (2) 

Magnetic    Resonance    Imaging 

The  phenomenon  of  nuclear  magnetism  was  well  known  long 
before  the  development  of  MRI  as  a  tool  for  the  examination  of 
morphological   structure   and   was   applied   to   the   more   basic   physical 
and  chemical  analysis  problems.     The  local  magnetic  field  of  a 
nucleus   is  responsible   for  the  exact   Larmor  frequency   of  the 
particular  spin  moment.     This   is  described  by  the  chemical   shift 
which   is   the   frequency   difference   of  a  particular  nucleus   compared 
to   some   standard.      This   property   along   with   many   other  more 


complicated   interactions   involving   steric   hindrance   and   thus 
relaxation  time  differences,   and  coupling  between  nuclei   through   a 
number  of  different  mechanisms,   allow   precise   characterization   of 
local  chemical  environments   and  thus   chemical   structure   (3). 

The  basis  of    MRI  is  to  eliminate  the  fine  spectral  details  and 
observe   only   the   nuclear   distribution   without   regard   to    small 
environmental   differences   between   nuclei.      The   common 
experiments   produce   contrast   through   differences   in   proton   density 
and  Tl   and  T2  differences,  although  sequences  for  many  other 
contrasts  have  been  developed  (4).     This  is  accomplished  in  the 
standard   imaging   sequences   by   producing   a   term   in   the   spin 
Hamiltonian   which   dominates   the  dynamic   behavior.     To   produce   a 
one-dimensional  picture   of  the   nuclear  density   in   a  sample,   a   strong 
linear  gradient  of  magnetic  field  is  added  to  the  homogeneous  main 
field.     This   serves  to  spread  the  frequency  to   such  an  extent  that 
chemical   shift  differences   are   normally   unnoticable.      Furthermore 
the   frequency   differences   encode   the   spatial   dimension   and   allow   an 
image  of  the  material  to  be  produced.     Application  of  a  radio 
frequency   (rf)  pulse   followed   by   acquisition  of  the   signal   induced   by 
the   moment  precession   during   application   of  a   linear   magnetic 
gradient  results   in  a  one-dimensional   image   which  is   a  projection  of 
the  sample  onto  one  line.     This  is  the  simplest  imaging  pulse 
sequence,   but  is   representative   of  one   dimension   in   almost  all   state- 
of-the-art    imaging    pulse    sequences. 

The   most  basic   imaging   sequence   currently   used   clinically 
involves   the   use   of  the   gradient   described   above   for   frequency 
encoding   in   one   dimension   and  then   a  repetitive   gradient   which 


varies  incrementally  with  an  acquisition  at  each  step.     This  gradient 
set  produces  phase  encoding   and   allows   an   inversion   via   Fourier 
transform  to  the  two  spatial  dimensions.     An  important  point  to  note 
for  future  reference  is  that  one  dimension  is  a  frequency  dimension, 
the  bandwidth   of  which   is   dictated  by   the   strength   of  the   frequency 
encode    gradient. 

The   NMR   Instrument 

The  introductions  to  NMR  and  MRI  above  refer  indirectly  to 
many  of  the  basic  sections  of  the  NMR  instrument.     The  primary 
requirement   is   a   magnetic   field.      Although   experiments   have   been 
carried  out  in  the  earth's  magnetic  field  (5)  and  fringe  fields  (6),  the 
standard    instrument   includes    a   carefully    designed   magnet   with    a 
strong,   homogeneous  magnetic   field.     The  magnet  may  be  a 
permanent    magnet,    an    actively    maintained    "resistive"    electro-magnet 
or  a  superconducting  electromagnet.      An  accessory  to  the   magnet  is 
the  shim  winding  set,  which  is  a  group  of  coils  which  can  be  driven  to 
correct   for  inhomogeneities   in   some  particular  region   of  the   magnetic 
field.   Closely   associated  are   gradient  windings   which  are  used  to 
produce  the  linear  magnetic  fields  used  in  imaging.     Typically,  the 
computer   which   controls   the   pulse   sequence,    sends   the   voltage   shape 
required   to   produce   the   current   (and   thus   gradient   field)   of  interest. 
These   three   components,   the   magnet,   the   shim   coils   and   the   gradient 
coils   along   with   their   associated   power   supplies   complete   the   static/ 
low    frequency    magnetic    field    requirements. 


The  next  basic  requirement  is  the  ability  to  perturb  the  spins 
from  equilibrium.     This  is  normally  implemented  with  the  use  of  a 
large  excitation  coil  tuned  to  resonance  at  the  Larmor  frequency  and 
a  high  power  rf  amplifier,  although  local  transceive  coils  are 
sometimes   employed  at  lower  power  levels.      For  modern   pulse 
sequences,   it  is  necessary  to  produce  shaped  rf  bursts.     As  with  the 
gradient   waveform,    the   computer   sends    information   to    a   digital-to- 
analog    converter   whose   voltage   modulates    the    Larmor   frequency 
sine  wave.     The  rf  coil  may  receive  as  well  as  transmit,  or  separate 
coils  may  be  employed.     In  either  case,  the  spin  magnetic  moment 
induces  a  voltage     in  the  coil  which  is  delivered  to  the  receiver. 

The.  NMR  Receiver 

Since  NMR  and  thus  MRI  inherently  generate  low  level  signals, 
it  is  crucial  to  produce  as  little  extra  noise  as  possible  during  the 
process   of  detection  and  recording  of  the   signal.     The   signal-to-noise 
ratio  (SNR)  is  the   standard  quantitative  measure  of  the  quality  of  a 
signal  in  the  presence  of  noise  and  will  be  discussed  below  and  in 
chapter  2.     The  block  diagram  of  the  standard  MRI  receiver  is  shown 
in  figure   1-2.     The  first  item  in  the  chain  is  the  probe  which  is 
sensitive  to  a  magnetic  field.     Typically  this  is  a  loop  antenna 
designed  to  be  sensitive  to  near  field  signals.     The  system  is  typical 
of  communications   receivers,   which   deal   with   an   audio   frequency 
signal   modulated   on   a   higher  frequency   carrier.      The   relaxation 
phenomenon    has    approximately    audio    frequencies,    while    the    Larmor 
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Figure  1-2.     Block  Diagram  of  NMR  receiver 


frequency  is  typically  tens  of  megahertz,  so  this  type  of  receiver  is 
well  suited  to  the  purpose.     It  is  important  to  note  that  the 
modulation   contains   both   amplitude   and   phase   information.   The 
incoming   signal   is   preamplified   and  eventually   mixed  down   to   audio 
for   digitization. 

The   signal-to-noise  ratio   is  the  measure  of  the  quality   of  some 
information  in  the  presence  of  noise.     Chapter  2  contains  specific 
discussions  of  formulations  of  SNR,  but  qualitatively     the  signal  is  the 
voltage   induced   in   a  receiving   prove  by   the   rotating   magnetic 
moment.     Similarly,  the  noise  is  a  random  voltage  induced  in  the 
same  probe.     There  are  many  possible  sources  of  noise,  but  in  the 
context   of  these   discussions  only   white   thermal   noise   will   typically 
be  considered.     This  is  related  to  any  dissipative  process  by  the 
fluctuation-dissipation   theorem   of  statistical   mechanics.      As   will   be 
discussed,   the   most   straightforward   quantitation   of  the   noise   is   in 


10 


terms  of  an  equivalent  electrical  resistance  of  the  probe  which  is  the 
addition  of  resistances  from  a  variety  of  sources. 

SNR  is  typically  measured  in  a  MRI  image  (which  is  a 
magnitude   Fourier  spectrum)   using   the   algorithm   depicted   in   figure 
1-3. 


Signal  and  noise 
mean  here 

Noise  mean  and 
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Figure   1-3     A  depiction  of  a  simple  image  of  a  round  phantom  with 
indications   for   SNR   measurement 

(sTn)-N 
The  SNR  for  the  image  depicted  is   SNR=       ■    ,     ,  „■     The  means  and 

std.dev(N) 

the   standard   deviation   are   with   respect  to   a  large   number  of  pixels. 

The   area  of  the   noise  which  is  used  for  the  measurement  is   arbitrary 

in  most  cases,   since  the  noise  is  uniformly  distributed  in  the  image 

space,  except  for  a  few  specific  cases  discussed  in  chapter  6. 

As  is  well  known,  noise  figure  (7)  is  a  standard  means  of 

characterizing   the   extra   noise   inserted   into   the   source   noise   and 

signal  by  system  electronics.     An  expression  for  noise  figure  for  a 

multi-stage    system    is 


G,  G,  G? 


1 1 


where  NF  stands  for  noise  figure  of  a  particular  stage  and  G  stands 
for  gain  of  that  stage.     It  is  clear  from  this  expression  that  the  first 
stage  is  critical  to  the  system  performance.     The  combination  of  low 
noise  figure  and  high  gain  in  the  preamplifier  is  most  desireable. 
Currently,   the   state-of-the   art  preamplifier  has   approximately   0.5   dB 
noise  figure  and  30  dB  gain. 

Figure   1-4  shows  the  typical  method  for  detection  of  the 
narrowband    signal   modulated    on    the    Larmor   frequency    carrier. 


<% 
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Figure  1-4.     Quadrature  phase  detection  receiver 


The    signal   is   first   down-converted   to   an   intermediate   frequency. 
This  prevents  the  existence  of  large  signals  at  the  Larmor  frequency 
which  could  affect  the  signal.     This  is  split  into  two  equal  signals. 
Each  signal  is  multiplied  by  a  reference  signal  at  approximately  the 
intermediate   frequency   with   the   phase   of  the   two   reference    signals 
maintained    90°   different.     The  low  frequency  component  of  the 
outputs   describe   a   vector  which   rotates   and   decays   with   time. 
Mathematically   only   a   single   down-conversion   is   equivalent   to   the 


1  2 


method  depicted  and  the  equations  for  direct  mixing  are  shown 
below.     A  particular  frequency  component  of  the  incoming  signal  is 
ce~"Tl  cos(coor  +  0).     The  process  of  quadrature  phase  detection 

produces 

Low  pass{ce~'/T'  cos(ay  +  (j))cos(corr)} 
and  Low  pass{ce~',T'  cos((u0f +  (t>)sin(corr)).    Using 

trigonometric   identities   these   can   be   rewritten   as 

Low  pass[ce~"T*  [cos((coo  +  cor)f  +  <)>)  +  cos((coo  -cor)r  +  <!>)]} 
and  Low  pass{ce"IT'  [sin(((0o  +  cor )r  +  <|>)  -  sin((coo  -  co, )t  +((>)]} . 
Typical   T2  decay  has  no  influence  on  the  separation  of  the  sum  and 

difference   frequencies   and  the   resultant   signals   are 

ce~'ln  cos((coo  -co,)f  +  <t>) 
and  -ce~'ITl  sin(((0o  -  CO,  )f  +  (j)).     If  these  signals  are 

assigned  to  the  x  and  y  axes  respectively,     a  spiral  pattern  is 
produced  which  traces  out  a  direct  representation  of  the  end  of  the 
magnetic   moment   vector   as   it   decays.      The   term   quadrature 
detection  refers  to  the  use  of  two  reference  signals  and  detectors. 
Single  detection  produces  only  a  projection  of  the  spiral  on  to  one 
axis.     This  results  in  ambiguity  as  to  whether  the  signal  is  greater  or 
less  than  the  reference  signal.     Besides  incomplete  information  of  the 
frequency  of  the  signal,  this  also  results  in  noise  which  is  both  above 
and   below   the   signal   added   together   in   the   resultant   acquisition. 
Typically  the  noise  in  NMR  is  due  to  the  source  probe/sample  system 
and  receiver  losses.     This  implies  that  the  noise  is  characterized  as 
white   noise   and   the   resultant   standard   deviation   is   thus    V2    higher 
than  if  quadrature  detection  is  utilized.     This  improvement  in  SNR  by 
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-s/2   due  to  the  use  of  quadrature  detection  is  of  interest  for  future 
discussion   of  quadrature  coil   technology. 

Original   Work 

The   above   descriptions   serve   as   basic   introductions   to   the 
concepts  to  follow.     Chapter  2  gives  more  detailed  background  on  the 
topics  necessary  to  the  original  analyses  in  the  remaining  chapters. 

The   bulk  of  this   dissertation,   chapters   3-8,   describes   the   results 
of  applied  physics   aimed   at  the   central   theme   of  improving 
sensitivity  in  NMR  and  MRI.    This  is  of  crucial  importance  in  the  field, 
because  the  source  signals  are  so  small  as  to  be  comparable  to  the 
thermal  noise  of  a  probe.     The  results  of  almost  every  interesting 
experiment  are  improved  by  an  increase  in  the  effective  SNR.     MRI 
continues  to  be  a  highly  signficant  clinical  tool  for  diagnosis  of 
disease  (8-10)   and  promises   new   utility  in  the  form  of  diagnosis  of 
breast  cancer  (11,12)  and     functional  imaging  of  the  brain  (13,14)   to 
name  a  few.     Even  when  increasing  SNR  does  not  reveal  new  detail  of 
structure   (morphological   or   spectral),   it   always   has   the   capacity   to 
improve    time    resolution. 

Specifically,   this  work  details  the  attempt  to  improve  the  SNR 
of  images  acquired  using  MRI  systems  in  the  range  of  0.5T  to  2.0T 
and  with   fields  of  view   approximately   5-40cm.      The   general 
philosophy    was   to   design   the   optimum   receiver   equipment   with 
knowledge  of  the  signals  to  be  obtained.     Though  similar  in  intent  to 
optimal   receivers   for   digital   systems,   the   information   available    about 
the  analog  signals  to  be  received  is  much  less  complete.     However, 
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some  improvements  were  made  over  the  state  of  the  art  for  general 
systems    and   further   improvements   were    made   for    specific    situations. 
The  system  from  the  coil  to  the  first  receiver  stage  were  considered 
for   improvement.      Sometimes   these   improvements   were   identified 
with  a  particular  type  of  experiment.     Again,  this  allowed  the   ability 
to   increase   specificity   and   thus   make   an   improvement   in   the   system. 

Because  this  work  has  been  in  an  extremely  active  research 
area,  some  of  the  topics  have  been  reported  in  the  literature 
following   the   authors   independent   developments.      Specifically 
sections  of  chapter  6  have  been  reported  by  other  workers  in  the 
field.     Unless  otherwise  noted,  it  can  be  assumed  that  all  of  chapters 
3-8   describe   the   independent   work   of  the   author,    with   invaluable 
assistance  and   support  from  the  UF  NMR  research  group. 


CHAPTER  2 
REVIEW  OF  STATE  OF  THE  ART 

Signal-to-Noise    Ratio 

Since  the  fundamental  theme  of  this  work  is  enhancement  of 
the  signal-to-noise  ratio,     a  model  is  needed  which  describes  the 
parameters  of  interest  in  SNR.     A  standard  formula  for  the  SNR  (1,15) 
following  a  90°  pulse  is 

SNR  =  K^°  Jllj '  [2J1 

K  is  a  numerical  factor  dependent  on  coil  geometry;  n  is  the  filling 
factor  (the  ratio  of  sample  volume  to  coil  volume);  M0  is  the  nuclear 
magnetization;    u0  is  the  permeabilty  of  free  space;  Q  is  the  quality 
factor  of  the  coil;  w0  is  the  Larmor  angular  frequency;  Vc  is  the 
volume  of  the  coil;  F  is  the  noise  figure  of  the  preamplifier;  k  is 
Boltzmann's    constant;    Tc  is  the  temperature  of  the  probe;  and  Af  is 
the  bandwidth   of  the  receiver. 

The  description  of  SNR  in  NMR  was  simplified  by  Hoult  and 
Richards  in   1976  (15).     The   improvements   over  the  older  formula, 
shown  above,   are  primarily  based  on  the  use   of  reciprocity  (a 
theorem   which   relates   transmission   performance   of  a   probe   to   its 
receive    performance)    and   the   identification   of  the    important 
variables   which   can   be   most  easily   controlled   or  at   least   measured. 
The   formula   is   presented   as 
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in  which  the  numerator  describes  the  signal  induced  by  the  spin  and 
the  denominator  describes   the   thermal   noise   generated  by   a   resistor 
(the  coil).     The  most  essential  part  of  the  equation  for  our  purposes  is 
the  Mo,  the  magnetic  moment,  is  equivalent  to  the  field  of  a  dipole 
rotating  in  a  plane  at  the  Larmor  frequency,  and  the  magnetic  field 
appears  in  the  numerator,  while  the  square  root  of  the  product  of  coil 
resistance  and  bandwidth  appear  in  the  denominator.      These  are  the 
only  parameters   that  will   be   used  in   the   following  discussions. 

In   1979  Hoult  and  Lauterbur  extended  their  discussion  of  SNR 
to  biological  tissue  (16).     The  expression  above  is  still  fundamental; 
however,  the  resistance  in  the  denominator  is  shown  to  come  from  a 
number  of  sources  of  relevance  in  MRI.     Reciprocity  is  here  applied 
to  NMR  probes;  however,  the  basic  concept  is  well  known  from 
general   electromagnetics.      The   primary   result  is   that   a  receiving 
probe   may   be   described   in   terms   of  its   transmission   properties. 
Using  reciprocity,  one  finds  that  the  power  dissipated  in  a  coil  can  be 
lost  in  4  fundamental  ways.     First,  the  coil  may  be  made  of  a  non- 
superconducting   material   which   implies   ohmic   losses   in   the 
conductor.      This   is  worsened  by  the   skin-effect  conduction  pattern   of 
current  at  radio  frequencies.      The   skin   depth   decreases   as   f_1/2 
(where   f  stands   for   frequency)   which   implies   that   the   resistance 
increases  as  f1/2.     Since  the  signal  induction  voltage  increases  as  f2 
(M0  increases   linearly   with   B0   and  thus  frequency,  and  the  time 
derivative    introduces    another    linear    dependence    of   frequency),    the 


SNR  increases  as  f14  for  ohmic  losses  only.  Second,  there  will  be  non- 
conservative  electric  fields  associated  with  the  coil  components.      If 
the  sample  is  conductive  these  electric  fields  will  deposit  power  in 
the   sample.      The   frequency   dependence   for  this   effect   is   complicated 
and  may  involve  a  resonance  (17).   Third,  the  magnetic  field  which  is 
the  primary  property  of  the  coil  will  also  generate  eddy  currents  in  a 
conducting  sample  which  represent  power  loss.     The  inductive  losses 
are   such  that  the  resistance  increases  as  f2  (16)  ,     thus  providing  a 
linear  increase  of  SNR  with  frequency  for  inductive  losses  alone. 
Fourth,   the  coil  has  a  certain  propensity   to  radiate  energy   as   an 
antenna.     This  is  directly  related  to  the  size  of  the  structure 
compared   to   a  wavelength   of  the   frequency   of  interest   and   the 
geometry   of  the   structure.      Radiation  resistance   of  a   simple   Hertzian 
loop  increases   as   f4  (18),  and  it  would  appear  that  there  is  no 
improvement  in   SNR  with  frequency;   however,   the  signal   induction 
formula  must  be  modified  as  well  if  the  radiation  is  an  important 
term. 

The   first  effect  can  be   addressed  by   maximizing   the  cross 
section  of  the  conductor  as  far  as  is  reasonable,  using  the  best 
available   material,   and   choosing   high   quality   discrete   components. 
For  mid  to  high  field  MRI,  this  loss  is  not  normally  dominant.     The 
second  effect  can  be  reduced  by  use  of  balanced  matching  (19), 
reducing  the  effective  inductance  of  the  coil   (20),   and  the  use  of 
Faraday   shields   (17).      All   of  these   methods   reduce   the  effective 
electric  field  in  the  lossy  sample.       An  example  of  the  process  is 
discussed   below   which   demonstrates   the   effect   of  a   particular 
implementation   of  reducing   the   inductive   reactance   of  the   coil. 


The  power  lost  due  to  the  electric  field  in  the  sample  is 
P=  [ITerjEfrfV,  where  a  is  the  conductivity  of  the  volume  and  E  is  the 


electric   field.     Figure   2-1    shows  the  voltages  on   the  components 

which    are   responsible   for   the   non-conservative   E-field    losses. 
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Figure  2-1.     A  comparison  of  two  tuned  circuits  with   different 
placement   of  components,   which   produce   different   power   losses 


Since   V=  [e  dx,  it  is  clear  that   1/4  as  much  power  is  lost  to  the 

sample   if  the   inductive   reactance   is   reduced   to    1/2   its   previous 
value. 

The  third  effect  remains  if  one  couples  to  the  spin  system.     The 
ideal  would  be  to  have  the  magnetic  field  only  in  the  region  of  the 
sample  of  interest  and  only   in   the   transverse  direction.      (The  rf 
magnetic  field  parallel  to  the  static  field  produces   loss  but  does  not 
couple    to    the    transverse    magnetization.) 

The   fourth  effect  can  be   moderated  by   use  of  geometries   which 
have    lower    radiation    resistances,    for   example    quadrupoles. 


1  9 


Distributing  capacitors  around  the  inductance  may  also  lower  the 
radiation  resistance  by  prohibiting  a  buildup  of  phase  of  current 
around   the   loop. 

There  is  one  potential  complication  in  the  equation  for  SNR 
given  above.     The  calculation  of  the  rf  magnetic  field  may  not  be 
accurate  if  done   in  the  absence  of  the  sample.     Normally,  however, 
the   permeability   of  biological   tissue   is   approximately   that   of  free 
space  and  the  conductivity  is  not  so  high  as  to  induce  currents  which 
are  on  the  order  of  the  coil  currents.     If  these  conditions  apply  and 
the  sample  dimensions  are  not  close  to  a  wavelength,  then  a  free 
space  calculation  of  magnetic  field  is  sufficiently  close.   (19) 
Radio-frequencv    Coils 

The  first  level  of  application  of  knowledge  of  the  problem  to 
improving  the  SNR  for  an  image  is  to  identify  a  typical  sample  and 
characterize  it  in  terms  of  size  and  electrical  properties.  For  optimal 
performance  a  receiving  probe  should  be  designed  so  as  to  pick  up 
only  noise  which  is  inseparable  from  the  signal  of  interest,  i.e. 
inductive  losses  only  in  the  region  of  interest  with  as  little  field  in 
the  static  field  direction  as  possible. 

The  first  NMR  based  systems  were  small  bore  magnets  for 
characterization   of  physical   properties   and   chemical   analysis.      The 
magnetic  probe  was  a    volume  coil  such  as  a  saddle  coil  (15)  or 
solenoid.   This    was   appropriate   to   match   the   typically   non-electrically 
conducting   small   samples   in   test   tubes. 

In  the  early  1980's  the  process  of  matching  the  coils  to  specific 
purposes  began  to  be  applied  to  clinical  imaging.  Surface  coils,  which 
were   basic   Hertzian    loops,   were   used   for  obtaining   inhomogeneous 
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but  high  SNR  images  of  structure  near  the  surface  of  the  human 
body.  The  first  use  of  a  surface  coil  in  vivo  was  described  by 
Ackerman  et  al  in   1980  (21).     This  work  dealt  specifically  with  in- 
vivo  spectroscopy;   however,   the  result  that  higher  local   SNR  was 
obtained   with   surface   coils   was   soon   applied   by   those   interested   in 
MRI.     Ackerman  et  al.   showed  the  well  known  magnetic  field 
produced  by  an  electrically  small  circular  loop  with  reference  to 
figure    2-2. 


Figure  2-2.     A  circular  current  loop  with  radius  a,   and  unit  current 


The  two  components  of  the   magnetic  field  are 

and 
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where  K  and  E  are  complete 


elliptic   integrals   of  the  first  and  second  kind. 
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By    1984   several  papers   dealing   with  region-specific   surface 
coils  had  reached  the  MRI  literature  (22,23,24).     These  papers 
quickly  showed  that  images  of  many  anatomical  regions  could  be 
greatly  improved  by  this  procedure.      From  the  description  of  the 
signal-to-noise  ratio  above,   it  was  clear  that  the  more   local   the  field 
to  the  region  of  interest,  the  lower  the  overall  resistance  of  the  coil 
and  therefore  the  better  the  SNR.     An  analytical  solution  for  the 
optimum  diameter  of  a  circular  surface  coil  was  described  in   1985 
(25),   and  was  later  taken   to  derive   the   ultimately  acheivable   SNR  for 
a  given  depth  (26). 

Before  this  time  the  only  coils  available  were  volume  coils  for 
the  whole  body  or  perhaps  the  head.     These  coils  include  the   saddle 
coil   (16)   and  various  versions  of  the  birdcage  resonator  (27).   The 
saddle  coil  is  derived  by  finding  the  positions  of  four  currents  lying 
on  a  cylinder  which  flow  parallel  to  its  axis  and  provide  the  best 
homogeneity  of  the  field  in  the  region  of  the  center  of  the  coil  (28). 
The  birdcage  coils  are  a  general  class  of  distributed  phase  coils  which 
are  analogous  to  a  full  wavelength  of  transmission  line  connected  in  a 
continuous  loop.     One  of  the  early  versions  of  this  was  exactly  this, 
where   the   transmission   line  was   coiled   into   a  toroid   shape   and   the 
center  conductor  was   exposed   on  each   loop  near  the   sample   (29). 
The   Alderman-Grant   resonator   is   a   two-leg   version   of  the   discrete 
element  implementation  (30).     The  general  version  of  this  coil 
described   first  by  Hayes  et  al.  (27),  may  have  any  number  of  legs 
and   the   current  in   these   legs   follows   a   sinusoidal   distribution   which 
has  been  shown  to  be  the  optimum  coil   for  cylindrical   geometry  at  a 
point   in   the   center  (31). 


About   this   time   quadrature   volume   coils   also   became   useful   for 
clinical  imaging  (32).    Quadrature  excitation  uses  half  the  power  for 
the  same  rotating  magnetic  field.       The  method  of  production  of  the 
quadrature   rotating   field   is   exactly   analogous   to   using   Euler's 
identity   to   produce   an   imaginary   exponential   with   the   sum   of  two 
sinusoids.      One   sinusoid   produces   a  rotation   and   a   counter-rotation 
each  with  an  amplitude  of  one  half.     Putting  all  input  power  into  the 
rotation  which  couples  to  the  spin  is  most  efficient  as  can  be 
observed  from  eq.   2.2.     The   arguments  demonstrate  lower  power  for 
excitation  and  reciprocally  higher  SNR  in  reception.     Thus  quadrature, 
like   the   quadrature   detection   receiver,   was   designed   to   receive   two 
equal   magnitude   signals   which   differed   in   phase   by    7t/2. 

The   birdcage   resonator   naturally   has   two   modes   related   to   the 
direction   a   traveling   wave   propagates   around   the   structure.      Driving 
the  coil  from  some  location  of  the  coil  will  non-preferentially  activate 
both  modes  with  the  result  that  the  two  traveling  waves  add  and 
form  a  standing  wave.     This  results  in   a  cosinusoidal  distribution  of 
currents.     At  a  position  90  degrees  from   the  drive  point  (where   a 
maximum    in    current    appears)    the    current   should    correspond   to 
cos(90°)  =  0.     Therefore  another  drive  point  can  be  located 
independently   at  this   location.     The   relative  phase   of  the   second 
drive  allows  cancellation  of  one  or  the  other  rotational  mode 
produced   by   the   first  drive   point   and   thus   will   determine   which 
rotational  direction  remains.      Any   other  volume  coil   is   made 
quadrature    simply    by    placing    another   coil    physically    perpendicular 
to  the  first  coil  and  driving  the  two  with  a  90°    phase    difference. 
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The  next  improvement  of  coils  was  the  application  of 
quadrature   principles   to  inhomogeneous  coils,   i.e.   surface   coils.     The 
first   reference   to   non   volume   quadrature   was   by    Arakawa   (33)   who 
described   the  principle   applied  to  slightly   curved   surface   coils.      The 
first  entirely  planar  quadrature   coil   was   described   by   Hyde   et   al.   in 
1987   (34).      Quadrature   surface   coils   have   since  become   quite 
common  with  routine  use  on  clinical   MRI  systems.     The 
generalization  of  quadrature  to  more  coils  and  other  phases   and 
magnitudes   is   discussed   in   chapter  5. 

Independent    Signal    Acquisition 

In    1986  Hyde  et  al.  (35)  first  published  a  discussion  of  the 
possibility  of  making  images  using  independent  coils.      In   this  paper  a 
claim   is   made   for  the   equivalence  of  independence   and   the   condition 
of  zero  mutual  inductance  for  coils.      Much  controversy  has   arisen 
concerning   this  contention.      Hyde  et  al.   presented   a   theorem   relating 
noise  correlation  and  mutual  inductance  in    1988   (36)  and  refined  in 
1991(37,  p.  36):   "If  two  coils  exhibit  zero  mutual  inductance,  there 
can  be  no  correlation  of  noise."     This  position  was  refuted  in   1989 
(38,  p.  402):  the  claim  that  "two  coils  with  no  mutual  inductance  will 
always   have   uncorrelated   noise. ...is   not   correct   in   most   instances"   and 
also  in    1990   (39,  p.   208):    "these   calculations. ...contradict   assertions 
that. ...coils   with   zero   mutual   inductance   have   no   correlated   noise." 
The   authors   described   the   noise   correlation   coefficient   as   an   electric 
coupling  coefficient.   In    1992   (40,  p.   85),   a   slightly  different  result 
was  given:   "if.. ..there  is  no  voltage  cross  talk  between  the  coils,   then 
noise,  either  from  coils  or  sample,  will  not  be  correlated."     Chapter  4 
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contains  a  concise  refutation  of  the  claims   summarized  above  from 
references  36,   37   and  40. 

The   fundamental  concept  is  to  receive   signals   from   separate 
and   somewhat   independent   coils   and   preserve   the   signals    and   noise 
to  allow  pixel-by-pixel   combination  of  the  signals.     The  degree  of 
independence   in   the   noise   determines   how   much   the   SNR   improves 
by  addition  of  the  two  signals. 

This  work  led  to  the  phased  array  which  was  described  by 
Roemer  et  al.  in   1990  (39).       This  is  a  general  discussion  of  multi- 
channel coil  systems.     Each  coil  has  its  own  receiver.     A  linear  array 
of  loops  is  a  simple  model  to  consider.     The  first  coil  receives  signal 
from  a  region  primarily  beneath  it.     The  second  coil  receives   some 
signal  from  the  region  of  the  first  coil,  but  mostly  from  sample 
immediately  adjacent  to  it.      If  one  considers  having   made   two  images 
separately  it  should   be  clear  that   the  optimum   image   of  the  two 
signals  will  have  SNR's  higher  than  either  of  the  single  images  in 
every   location.      The   advantage   this   brings   over  quadrature   is   the 
ability  to  arrange  the  coils  in  the  z  direction  along  the  main  field  and 
to   combine   each   pixel   with   different   combination   weighting   to   obtain 
optimum   performance   everywhere.      It   should   be   pointed    out   that 
coils  which  are  coupled  to  the  same  lossy  sample  will   not  generally 
be  entirely  independent.     This   is   because   noise   generated   in   the   lossy 
medium  can  couple  to  both   coils  producing  a  correlation   which 
prevents    complete    independence. 

One  of  the  problems  solved  by  Roemer  et  al.  (39),  is  the 
apparent   necessity   to  decouple   the   receiver  coils   as   much   as   possible. 
Voltage   crosstalk   between   coils   introduces   a   noise   correlation   and 
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non-optimum    signal   combination   (in   general).      The    method    employed 
is  to  mismatch  the  coil  to  the  preamplifier  in  such  a  way  that  the 
current   induced   in   each   loop  during  reception   is  much   lower  than   it 
would  be     in  the  matched  condition.     In  a  signal-to-noise  sense,  this 
process   is   non-lossy.     The  reduction   in  current  means  that  the 
receiving  coil,  as  a  source,  loses  efficiency,  resulting  in  less  transfer 
of  signal  and  noise  to  another  receiving  coil.     The  intrinsic  noise 
correlation  induced  by  the  sample  is  not  influenced  by  this  method. 
It  is  shown  in  chapter  4  that  isolation  is  not  actually  necessary  for 
optimization  as   long   as   the  coupling   is   measured  under  experimental 
conditions. 

Another  topic  discussed  in  the  work  of  Roemer  et  al.   (39)  is  the 
problem   of  image   reconstruction   utilizing    the   multiple   images 
received.      The   general   conclusion  is  that  the   optimum   treatment  of 
the   images   requires   knowledge   of  field   characteristics   of  the   coils 
used  as  receivers   and  noise  correlations   between   each  pair.      A 
complex   addition   utilizing   this   information   gives   the   optimum   SNR 
for  each  pixel  in  the  final  image.     As  a  matter  of  practical  importance, 
a  much  simpler  algorithm,  i.e.  the  square  root  of  the  sum  of  the 
squares   of  the   individual   image   pixels,   give   a   nearly   optimum   result- 
-   at  worst  about   10%  inferior. 

There   are   many   ways   in   which   any   intrinsic   independence   is 
preserved  while  still  getting  all  of  the  signals  to  the  computer.     The 
most    general    is    described    above— using   independent   receivers    all    the 
way   to   the   a-to-d   converters.      Another  method   which   will   be 
described  in  chapter  6  is  by  time  multiplexing  the  data  (41).     The 
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requirements   for  this   will   be   discussed   later  but   the   result   is   exactly 
the  same  as  for  the  multiple  receivers. 

Another   method   is   a  more   specific   procedure   in   which   narrow 
band   filters   at   the   Larmor  frequency   are   matched  to   the   coil 
positions   and  field  characteristics   (42).     The  great  benefit  of  this 
method  is  that  the  output  can  be  fed   to  a   standard   single  channel 
system   and  no  software  modifications  are   necessary.        This  will   be 
further   discussed   in    chapter   6. 


CHAPTER  3 

VARIABLE  FIELD  OF  VIEW  COILS 

Introduction 

A  variable  field  of  view  surface  coil  was  designed  to  allow 
variation  of  the  field  of  view  to  be  imaged  while   maintaining  the 
tuning     of  the  coil  to  a  substantial  degree.     This  is  an  attempt  to 
improve  imaging  by  allowing  a  significant  size  region  to  be  imaged 
but  still  enable  the  user  to  obtain  optimal  performance  on  a  smaller 
region.     This  work  was   completed  prior  to  the  appearance  of  multi- 
channel   systems. 

Radio  frequency  (rf)   surface  coils   have  become  a  necessity  for 
high  resolution  MRI  of  the  human  body.     This  is  because  optimization 
of  the  SNR  primarily  requires  that  the  size  and  shape  of  the  receiver 
coil  provide  an  optimum  filling  factor  given  the  region  of  interest 
(43).     Small  rf  coils  provide  larger  Bi   fields  per  unit  current  while 
presenting  a  smaller  surface  area  to  the  load.     The  reduction  in  the 
surface  area  of  the  load  results  in  less  degradation  of  the  coil  circuit 
Q,  reducing  coil  losses  (17,44,45). 

Popular  regions  of  interest  for  surface  coil   designs   include   the 
orbits,  TMJ,  neck,  spine,  heart  and  extremities.     A  number  of 
different    geometries    with    varying    degrees    of   performance    (46-50) 
have  been  proposed  to  optimize  SNR  for  each  of  these  regions  of 
interest.     Due  to  a  large  demand  for  MRI  of  the  nervous  system,  this 
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work  has  focused  on  the  spine.     This  area  is  especially  challenging 
because   the  clinician  would  like  to  survey  the  entire   spine  in   a   single 
acquisition  and  at  the   same  time  obtain  the  highest  possible   SNR  for 
high  resolution  from  thin  slices.     The  result  is  that  surface  coils  for 
the  spine  may  be   a  compromise   between   field  of  view   and   high 
performance. 

Several   workers   have   attempted   to    solve   this   problem   by 
constructing  coils   which  either  can   be  repositioned   (51)   or  have   sets 
of  elements  which  can  be  switched  into  place  to  change  the  field  of 
view   (52-54).      Repositioning   devices   can   be   very  effective   and  easy 
to  use  because  there  is  no  need  to  reposition  the  patient:   the  coil 
simply   moves   in   a   space   created   between   the   patient  and   the   table. 
However,   larger  fields   of  view  require   multiple  acquisitions.     The  use 
of  multi-element  coils,   however,   is   limited   by   the   number  of  sections 
which  can  be  effectively   switched.      In   addition,   such   switching 
requires   PIN   diodes   at  each  junction   which   may   reduce   the  efficiency 
of  the  design  by  reducing  coil  Q.     Another  approach  to  the  problem  is 
to  construct  an  array  of  coils  which  have  no  mutual  inductance  and 
to  connect  each  coil  to  a   separate  preamplifier,   phase  detector  and   an 
analog-to-digital  converter  (39).      This  approach  is  effective;   however, 
it   requires   much   additional   electronic   modification   of  the   receiver 
system   which   is   quite   expensive. 

The  objective  of  this  work  was  to  provide  a  flexible,  efficient, 
low-cost  solution  to  the  problem  of  varying  the  field  of  view.     The 
general  strategy  was  to  design  NMR  coils  that  would  be  variable  in 
size  but  remain  electrically   stable   during   adjustment.      This   kind  of 
probe   would   be   applicable   to  clinical   imaging,   where   tuning   and 
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matching    for   individual   patients   is    highly   inefficient   and    undesirable. 
A   common   method   of  adjusting   dimensions   of  electrical   systems   is 
the  use  of  the  trombone  (55).     This  adjustment  of  size  coincides, 
however,   with  a  drastic  change  in   the  impedance  of  the   system.      In 
fact,   trombones   are   typically   used   to   tune   to   a   particular   impedance 
and  have  even  been  utilized  for  this  purpose  in  NMR  (56).  The  goal 
was  to  make  a  self  compensating  trombone  which  adjusts  with  no  net 
change  of  impedance   (tuning  in  the  NMR  probe). 


Methods 
A  prototype  flat  rectangular  coil  was  constructed  that  could  be 
adjusted  from  about  12.5  cm  x   17.5  cm  to  about   12.5  cm  x  32.5  cm 
and  was  made  from   1/2"  and  3/8"   tubular  copper.     This  coil  is   shown 

in   figure   3-1. 

Inner 
copper  tubing 
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Dielectric  on 
middle  tubing 
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Outer 
copper  tubing 


Figure   3-1.      A   schematic   representation   of  the   adjustable   coil 
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Using  this  coil  as  a  model,  an  analysis  can  be  performed  which  will 
lead   to   an    approximately    constant   resonant   frequency    over   a 
considerable  range  of  size  adjustment.     This  analysis  assumes  that 
coil  dimensions  are  small  enough  to  ignore  wavelength  effects  and 
negligible  resistance  in  the  coil.     Figure  3-2  shows  a  circuit  model  of 
the  coil  and  figure  3-3   shows  a  simplified  version. 

Li        C(x)     L2+L(x) 
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Figure  3-2.     Circuit  model  of  the  adjustable  coil 
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Figure  3-3.     Simplified  circuit  model  of  the  adjustable  coil 


For  constant  tuning,    L      must  be  a  constant,   where 


;ff>Leq=;cfjLT(A')-;- 
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and  00  =  ■ 


For  the  case  of     a  simple  trombone,   schematically  shown  in 

figure   3-4,    the   following   equations   approximately   apply    (57): 
L(jc)  =  1.97E-  7 x[ln(4x/62 )  +  d,  fix  -  0. 75]  H 
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C(°-n^75^F 

where  linear  dimensions   are   in   meters   and   SI  units  are  used   for   £0, 
the   permittivity   of  free   space   and    £r,    the   relative   permittivity    for   the 
dielectric. 
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Figure  3-4.     A  schematic  of  a  trombone  section 


Also  CT(x)  =         1    (X)       and  LT(x)  =  2[L  +L,  +L(x)l. 
C(jc)  +  2C,  T  L   1       2  J 

In  order  to  satisfy  the  constraint  that   Leq  be  constant,  C(x)    and/or 

L(x)  must  be  variable.     As  can  be  seen  from  the  equations  above,  if 
the  trombone  is  used,  there  is  only  one  way  to  make   Leq  a  constant. 

This  is  to  make  er  vary  as  a  function  of  x  to  give  LT(x)   and 
l/(X)2CT(x)   equal   derivatives   in   x.     By  slightly  modifying  the 
trombone,   the  goal   of  constant   Lsq  can  be   attained  in   several 

different    ways. 

Results 
The   approach   presented   in   the   following   section   leaves   the 
inductive   change   nearly   the   same   as   in   a   standard   trombone   to 
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simplify   analysis.      For  the  protoype   discussed,   the  capacitance   of  the 
system    decreases    too   rapidly    (the    resonant    frequency    increases)    as 
the  coil  is  enlarged.     Two  methods  were  utilized  to  compensate  for 
this    situation. 

The  first  method  utilizes   a  set  of  telescoping   conductors   similar 
to  the  radio  antenna  of  an  automobile.     From  figure  3-5,  it  can  be 
observed  that   the  rate   of  decrease   in   capacitance   is   slowed   with   this 
method. 


Figure  3-5.     A  set  of  telescoping  conductors 

One  can  vary  the  diameter  and  length  of  each  segment  to  make  the 
resonant    frequency    approximately    constant    (within    the    Q    bandwidth) 
over  the   range   of  size   variation   desired. 

The  second  method,  shown  in  fig.  3-6,     involves  the  use  of  a 
layer  of  copper  foil  between  layers  of  dielectric.     This  foil  can  be 
easily   trimmed   to   compensate   for  the   fast   decrease   in   capacitance. 

The   two   methods   described   can   be   used   separately   or 
combined   to  obtain   the   desired  result.      Figure   3-7   shows   empirical 
results   for  the  prototype   coil   described   using   a   normal   trombone,   one 
sliding  inner  conductor,  and  an  inner  conductor  plus  a  piece  of  foil 
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Figure   3-6.      A  trombone  modified  with  extra  layers  of  conductor  and 
dielectric 


tape  of  a  particular  shape.     The  results  demonstrate  that  the  change 
in  resonance  frequency  with  change  in  coil   size  can  be  made 
negligible   with  appropriate  modifications.      Similar  results  could  be 
obtained  for  coils  made  with  foil   strips  for  the  conducting   surfaces. 

To  evaluate  the  coil's  effectiveness  in  a  clinical  setting,  a 
volunteer  allowed  imaging  of  his   spine.     It  is  frequently  of  interest  to 
obtain  a  large  field  of  view  of  the  spine  and  then  to  concentrate  on 
one  or  two   vertebral   bodies   in   particular.      Figures   3-8a   and   3-8b 
show  the  results  of  using  the  coil  in  an  extended  and  a  fully 
contracted   position   respectively.      The   images   were   obtained   at    1.5T 
with  a  GE  Signa  using  a  spin  echo  sequence  with  256  x   128  pixel 
resolution,   500  ms  repetition  time,   30  ms  echo   time,   four  averages,     5 
mm  thick  slices  and  a  40  cm  field  of  view. 
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Figure   3-7.      Empirical   results  of  resonant  frequency   vs.   coil  extension 


Conclusion 
The  possible  applications  of  the  basic  idea  are  numerous.     Any 
probe  which  has   a  linear  or  circular  segment  can  be  made  in  this 
way.     Also,  if  one  does  not  wish  to  attempt  to  modify  the  trombone,  a 
varactor  diode   can   be   used   to   maintain   a   constant  resonant 
frequency.     The   size   could  also  be   varied  remotely  with  hydraulic  or 
pneumatic    devices    or   motors. 
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Figure  3-8.  Spine  coil  images  a.  coil 
extended  b.  coil  contracted 
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Another   slightly   different   application    involves    any   case   in 
which  two  loops  are  coupled,  producing  splitting  of  the  resonant  peak 
dependent  upon  the  separation  of  the  two  loops.     By  utilizing  the 
above   modified   trombone   approach,   it   should   be   possible   to   vary   the 
coupling  while  forcing   the  peak  of  interest   to  remain   at   constant 
frequency.     An  example  of  this  is  an  adjustable  cardiac  coil. 

It   should   be   noted   that   the   resonant   frequency   does   depend   on 
coil  resistance  as  well  as  reactance.     If  the  Q  is  fairly  low  this  may  not 
be   negligible.      Therefore,   both   the   resonant   frequency   and   match 
cannot   be   exactly   corrected   because  the  extension   makes   the   coil 
have    greater    resistance. 

The   work  presented  here   was   previously  published  in   Magnetic 
Resonance   in    Medicine    13,   378-384  (1990)  and   forms  the  basis  for 
US  patent  #  5,049,821    assigned   to  the  University  of  Florida. 


CHAPTER  4 

COUPLING  BETWEEN  NMR  COILS 

Introduction 

Quadrature   combination,    as   introduced   in   chapter   2,    was 
implemented   primarily  for  use  with   volume   coils   and   provides   a   root 
two  improvement   in   SNR   and  half  the  power  required   for  excitation. 
Hyde   (34),   and   Arakawa   (33)   initiated   the   use   of  non-volumetric 
coils    with    quadrature    combination    and    substantial    improvements    can 
be  produced.     The  real  difference  is  that  the  magnetic  fields  of  these 
coils   are   not  homogeneous,   which   implies   that  the   optimal 
combination   occurs   only   for  those   regions   where   the   magnetic   fields 
of  the  two  coils  to  be  used  in  the  combination  are  of  equal 
magnitudes    and    are    perpendicular. 

The  theory  of  quadrature  is  insufficient  to  describe  all 
situations   of  different   field   magnitudes   and   arbitrary   crossing   angles. 
The  following  two  chapters  develop  a  theory  to  cover  all  situations. 
This   will   be   denoted   generalized   quadrature   and   is   also   accurately 
described  as  optimal   signal  combination  which,  however,   also 
includes  consideration  of  coupling.     The   general  topic  of  optimal 
combination  will  be  broken  into  two  sections.     In   this   chapter 
compensation  of  signal  and  noise  coupling  will  be  considered.  Chapter 
5  deals  primarily  with  the     generalization  of     phasing  to  the  array 
with    inter-coil    noise    correlation. 
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Since   the   introduction   of  the     quadrature  coil   (32),   the   first 
receiver  coil  array  in  MRI,  the  issue  of  isolation  (32,  58)  has 
frequently   been   an   area  of  concern   and  controversy.      Because   of 
existing   disagreement,   especially   on   the   relationship   of     noise 
correlation   and   coupling   (37,38,40)   as   discussed   in   chapter   2,   the 
general     problem  will   be  examined  by  reducing  it  to  more  specific 
areas  of  interest.     A  semantic  description  of  coupling  and  noise 
correlation  needs  to  be  carefully  given. 

cross    talk-this   can  also   be   referred   to   as   coupling-cross   talk 
between  two  coils  is  non-zero  if  and  only  if  a  current  originating  in 
coil  1   induces  a  current  in  coil  2. 

noise    correlation    (nc)-is   zero   if  and   only   if  the   noise 
voltages   delivered   to  each   coil's   load   (preamplifier)   can   be   added 
with   arbitrary   phase   with   the   result  that  the   total   noise  power  is   the 
same  as  the   sum  of  the   individual   noise  powers. 

intrinsic    nc  --noise    voltage    originating    from    a    shared    lossy 
load.     This  is  independent  of  crosstalk. 

extrinsic    nc  --noise    voltage    originating    from    the    coupling    of 
one  coil's  noise  voltage  to  the  other  coil 

It  is  shown  here  that  cross-talk  and  noise  correlation  are,  in 
principle,   separable.     Consider  the  circuit  model  in  figure  4-1.     The 
coil  matching  is  done  as  described  in  the  work  of  Roemer  et  al.  (39). 
The  coils  are  positioned  so  as  to  have  zero  mutual  inductance  in  the 
absence  of  a  load  and  it  is  assumed  that  the  load  does  not 
substantially   change   the   field   pattern.      Therefore   the   coupling   path 
between  coils  is  confined  to  the  path  through  the  load.     The  model 
shown   includes   each   coil's   isolated   noise   sources   into   the   resistance 
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Figure  4-1.     Equivalent  circuit  for  two  coils  loaded  by  the  same 
lossy   sample  and  thus  exhibiting  intrinsic   noise  correlation 
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in  series  with  the  coil  inductance.     Since  each  coil  is  coupled  equally 
to  the  model's   inductive   load,   this  resistance  produces  purely 
correlated    noise. 

To  examine  the  process  of  cross-talk,     assume  a  voltage   source 
va  in  series  with  the  inductance  of  coil   1   (vb   shorted  for  this  case). 
It  is  of  interest  to  know  how  much  the  induction  of  signal  and  noise 
into  coil   1   is  transferred  to  coil  2.     The  cross-talk  constant  will  be 
defined   as   i2/il-     It  is  easily  shown  that    lim  —  =  0  This  was 

demonstrated  by  Roemer  (39)  for  the  case  of  no  intrinsic  noise 
correlation. 

Now  consider  the  noise  which  is  generated  by  the  resistor  of 
the   lossy  inductive   load   (vb   active,  va  shorted).     It  is  clear  without 
any  calculation  that  the  induction  of  the  noise  voltages  into  coil   1    and 
coil  2  is  identical  and  not  a  function  of  Ramp.       The  noise  delivered  to 
amplifier    1    comes   from  the  inductive   load  resistor  and  the   local 
resistor  of  coil   1.  The   relationship   between   these   voltages   is   also   not 
a  function  of  Ramp-     The  conclusion  is  that  intrinsic  noise  correlation 
and    cross-talk    are    entirely    separable    whereas    cross-talk    is    entirely 
responsible  for  extrinsic   noise  correlation.      Any   lossless   method   of 
increasing   isolation   should   result  in   reduction   of  extrinsic   noise 
correlation.     Intrinsic  noise  correlation  is  impossible  to  change   in   a 
lossless  way.     With  these  factors  in  mind  it  is  logical  and  useful  to 
treat   noise   correlation   and   coupling   as   distinct   phenomena. 

The   remainder   of  this   chapter   discusses   the   investigation   of 
coupling   between   MRI  coils   and  the  effect  on   the   attainable   signal- 


to-noise  ratio  (SNR).         It  will  be  shown     that     if  the  signals  are 
combined  properly,  there  is  no  loss  in  SNR  due  to  coupling.     Long 
wavelength   approximations   will   be   used   and      it  will   be   assumed   that 

Isolation  (dBl 

the  degree  of  effective  coupling   (k  =   10        x      )  is  not  equal  to     one. 
In   chapter  5   the  issue  of  noise   correlation  and  the  generalization  of 
quadrature   gain   will   be  examined.      The   technique   of  reciprocity   (15) 
will  be  used  to  examine  SNR  by  assuming  the  coils  are  driven,     and 
the   principle  will   be   confirmed   with   a  comparison   to   signal   reception. 
Coupling   with   no   shared  resistance   will   first   be   examined   and   several 
examples   will   be   given   before   describing   the   general   situation. 

Discussion 
Relative  SNR  is  the  parameter  of  interest  and  this  is 
represented   via   reciprocity   as   the   rotating    component   of  magnetic 
field  divided  by  the  square  root  of     power  input  to  the  system.     This 
method  is  clarified  by  example.     Consider  the  single  saddle  type  coil 
with   arbitrary   arc   length   shown   in   figure   4-2. 
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Figure  4-2.     A  saddle  type  coil  seen  from   an  end  perspective 

A  point  of  interest  (POI)  is  defined  at  the  center  of  the  coil  .     The  coil 
has  some  resistance  R  which  in  general,  is  due  to  coil  electronics  and 
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local  sample  loading.     When  the  coil  is  driven  with  power  P,  the 
current  in  the  loop  must  be  of  magnitude   ^/P/R.     It  is  assumed  that 
any  matching  circuit  loss  can  be  lumped  into  this  loop  resistance. 
The  magnetic  field  is  a  linear  function  of  current  and  we  define   B  as 
the  amplitude  of  the  rotating  field  produced  by  unit  current     at  our 
POI.    B/VP    is  independent  of  current  for  a  single  coil  and  is  a  point- 
wise  measure  of  efficiency  exactly  analogous  to  SNR,  thus 

IrI         B;  R 

lBl  ■■-         I    --g-       [4.1] 

Vp    VtfR    ^R 

The  next  step  is  to  consider     two  coils  which  by  assumption 
have  no  intrinsic   noise  correlation.  This  implies  that  whatever  the 
mechanism   of  coupling,   there   is   no   resistor  other   than   the   individual 
isolated  coil  resistances.     To  see  the  effects  of  coupling,     the  case  of 
completely   independent  coils   are   first  considered.      Figure   4-3   depicts 
two   identical   reduced   arc   saddle  coils   which   have   been   oriented   for 
zero   mutual   inductance. 


Figure  4-3.     Two  reduced  arc  saddle  coils,  oriented  for  zero  mutual 
inductance 
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There     is  indeed  a  rotational  position  which  corresponds  to  zero 
mutual  inductance  aside  from  the  usual  position  of  a   7l/2  rotation.     To 
observe  that  this  is  so,  assume  that  the  coils  are  positioned  directly 
on  top  of  one  another.     They  will  be  strongly  coupled  with  all  of  the 
flux  from  one  coil  passing  with  the  same  sense  through  the  other.     If 
one  is  rotated  until  it  does  not  overlap  at  all,  the  flux  of  the  first  coil 
now  passes  through  the  other  loop  with  the  opposite  sense.     It  is 
expected,  in  practice,   that  the  flux  linkage  would  be  a  continuous 
function  of  angle  which  must  have  a  zero  between  a  positive  and  a 
negative  region.     Therefore,  it  is  possible  to  overlap  the  coils  for  zero 
mutual  inductance  and  the  angle  §  is  defined  as  the  physical  angle 
(between  bisection  axes  of  the  coils)  at  which  this  occurs.     This  angle 
also  gives  the  angle  at  which  the  field  lines  at  the  POI  cross. 

Assume  the  loops  have  the  same  resistance  R  and  coil  1  is 
driven   with   current    i,  and  coil  2  with  i2.    The  fields  cross  at  angle  § 

and     correspondence  of     the  x-  axis  with  real  numbers  and  the  y-axis 
with   imaginary   values      gives  the   total  rotating   magnetic   field 
efficiency: 


Bi,  +Bi2eJ 


B      i,  +  i2e 


i.|2R  +  kl2R      VR 


'11    -"    1*21    «  M     +\h\ 


[4.2] 


Since  the  fields  are  rotating  we  note  that  a  positional  delay  is 
equivalent  to   a  time   delay.      Suppose   the   currents   are  electrically 
varied  in  an  equal  power  split  so  that  one  of  the  currents  is  delayed 
with  respect  to  the  other  by  the  angle  0.     The  equal  power  split  will 
result   in   equal   magnitude   currents   and   thus   the   efficiency 
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B    1  +  eV 


It  is  clear  that  the  choice  of  6  = 


produces  the 


a/R       a/2 

result  that  the  total  rotating  magnetic  field  efficiency  is   a/2    times 
greater  than  for  the  linear  coil.     This  action  simply  time  delays  the 
field  which  is  positionally  advanced   so  that  the  fields  magnitudes 
add.   This   is   the   generalization   of  quadrature   gain  (59)   to   any 
physical   angle   given   the   important   assumption   that   the   coils   are 
truly  independent.     This  concept  is  more  fully  developed  in  chapter 
5.     It  should  be   noted  that  this  principle  directly  applies  for  low 
frequency  coils,  micro  imaging  coils,     coils  with   nonconducting 
samples   and   any   other   application   where   the   resistance   of  the   coil 
conductor  is   the  dominant   source   of  resistance,   and  there  is   no 
coupling   between   the   coils. 

Since  the  discussion  above  shows   that   a/2    gain  can  be  obtained 
for  any  physical   angle  if  the  coils  are  independent,  the  case  of 
coupled   circuits   is   next  examined.      The   description  will  remain 
general  by  assuming  the  change  due  to  coupling  can  be  described  as 
shown  in  Table  4-1. 


Table  4-1.     The  forms  of  coupled  vs.   uncoupled  coils  currents 


coil  1 


coil  2 


uncoupled 
coupled 


H 


I]  =  i]    +  ki2eJP 


12  =  i2    +  ki]eJP 
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The   parameter   p    describes   the  phase   angle   with   which   the   coupling 
occurs,  k  is  the  magnitude  of  the  coupling  and  6,  as  before,  is  the 
electrical  phase  angle  imposed  on  one  of  the  coils  by  the 
splitter/combiner.      Note   that   the   currents   which   are   seen   as 
originating  with  a  given  coil  are  not  the  same  as  the  uncoupled  coil 
currents  but  they  converge  for  small  k.     In  the  coupled  case  ij   is  the 
current  produced  in  coil  one  when  driven  and  with  the  other  coil 
present  but  not  driven.      Figure  4-4   shows   an  electrical   block  diagram 

of  the  coils  in  figure  4-3.       The  total  efficiency  of  the  pair  of  coils  with 

|BI,+BI,e*| 
the  net  currents   I,  and  I,  after  coupling  is  .        [4.3] 

t/|I,|2R  +  |I2|2R 

If  these   net  currents   can   be   manipulated   as   freely   as   the   uncoupled 
currents   by   variation   of  the   splitter  circuits   the   condition    I2  =  \^e~* 
may  be  satisfied  and     the  maximum  gain  remains    V2.     Note  that 

complete  freedom  of  manipulation  is  impossible  if  k  =    1.     The  splitter 

i  1-  ke^e* 

solution  is  given  by         —  =  e  '* — -.  [4.4] 

i,  1-  ke^e'^ 

This   implies   generally   different  magnitudes   of  the   driving   currents 
for  equal   amplitude   net  currents   after  coupling.      Note   that  the 
imbalance  in     /,  and  (',    does   not   correspond   to   amplification   but 
rather  a  shifting  of  power  from  one  loop  to  the  other.     No  active 
devices   are  required.      An  example  of  this  process  is  impedance 
matching  to  some  arbitrary  impedance.     This  scales  the  noise  and 
signal  equally  if  the  circuit  is  lossless. 

Because  the  ability  to  obtain  lossless  performance  in  the 
presence  of  arbitrary  coupling  may  be  unexpected,  the  gain  in 
efficiency   with   coupling  will   be  examined  in  detail   for  several 
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Figure  4-4.     Electrical  block  diagram  of  the  coils  of  figure  4-3 
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specific  cases.  The  single  coil  field  value,  B/VR,  will  be  separated 
and  the  remainder  designated  as  the  gain  G  to  which  we  will  refer 
hereafter.     The  total  gain  in  efficiency  is  then 


ll  +  AkeJI9+P)  +  (Aeje  +  keJ|i)e'*| 
G  = 


[4.5] 


||1  +  Ake'^f +\{Ae">  +  ke*)e»\ 
where    i2  =  Aefiit . 

Case   I   deals   with   parallel   fields   and   case   II   addresses   perpendicular 
fields. 


ll  +  Akey,e+P)  +  (Aefi  +  kejP) 
I.  0  =  0  G(<|)  =  0)  =    '  ,  '  [4-6] 

'    +  Ake'^H +\Ae"'  +  ke'p 


Selection  of  e  =  0  and  A  =  1  results  in  gain 


ll  +  ke;p  +  (1  +  ke®)\  ,- 

G(6f  =  0,9  =  0)  =    '  ,  ;    =V2.    [4.7] 

1  +  ke'p|  +|l  +  ke^l 


Note  that  this  particular  condition  leads  to  equal  currents  in 
the  loops  which  add  in     phase.       Figure  4-5  shows  vectors 
representing   the   rotating   fields   (at   some   instant   in   time) 
from  the  two  coils  at  the  POI  which  are  added  under  the 
conditions   described   above.   The   same   arguments    apply   for 
6  =  Jt. 


<()  =  0 

Independent  loops         Coupling,  arbitrary  (3      coupling,  9=0 


^^T     j1 j< 


Figure  4-5.      Vectors   representing   the   additions   of  the   rotating 
fields  described   in   case   I 

II.   $  =  71/2 

A.  Consider  the  special  case  p  =  0  with  A  =  1: 

ll  +  ke*  +  (e*  +  k)e""2\ 

G(<(>=  7t/2,p  =  0)  =         I '.     [4.8] 

V2{l  +  2kcos0  +  k2} 

If  one  chooses  the  splitter  angle  to  be  9  =  -7T./2,  the 

resultant  gain  is   ,1 .     However,  the  currents  in  loop  one 

Vl  +  k2  V 

and  loop  two  after  coupling  are  of  equal  magnitude  but  are 

not  out  of  phase  by  Jt/2.    The  angle  9    should  be  chosen  with 

full   knowledge  of  the  problem,   including   the   description  of 

coupling.     It  is  easily  shown  that  the  maximum  gain  of 

V2  occurs  at  the  angle  where  cos6    =   -    .    Fie.  4-6a 

1  +  k2 

graphically    demonstrates    the    field    addition    under    these 

circumstances. 


4') 


Independent  loops  coupling,  p=0  coupling,  6=-jt/2 


"x 


*i 


,    -2k 

y  =  cos 


\  +  k: 


N 


Figure  4-6a.     Field  additions  for  Case  II,  part  A 

It  is  important  to  note  that  the  splitter  angle  is  simply  the 
angle  which  produces  a  net  current  phase  difference  of       -n/2 

when  coupling  is   considered.      Again   the   same   argument 
applies  for  p  =  n. 

B.  Consider  the  special  case  p  =  n/2: 


For  A  =  1  and  6  =  -n/2,    a  gain  of  J is  obtained.     This 

'  6  Vl  +  k2 

result   is   the   same   as   above   and   has   been   previously   obtained 
for  standard  quadrature  coils   (2).     In  this  case  there  is  no 
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phase  which   produces   the   full   gain.     Another  degree   of 

freedom   must  be  exploited  -the   combiner  power  split.        On 

r             a        J  +  k 
optimization  one  finds  that  the  gain  is       V2  when  A  = — 

.     It  is  no  coincidence  that  this  equalizes  the  currents  in  the 
loops.     Figure  4-6b     shows  the  field  additions  for  these  cases. 

(j)=7t/2 

Independent  loops  coupling,  p=jt/2  coupling,  6=-Jt/2 

— ►/    4 —    — ►-►  — ^ 

v. 


coupling,  A  =  


Figure  4-6b.     Field  additions  for  Case  II,  part  B 

Reciprocity   predicts   the   above   results,   but   care   must   be   taken 
in  proper  consideration  of  the  situation.     The  special  cases  above  are 
illustrations  that     if  one  drives  two  equivalent  loops,   such   as   shown 
in   figure   4-3,   which   have   only   internal   resistance   (no   shared 
resistance)   and  there   is   no   other  power   sink,   then   forcing   equal 
magnitude   currents   with   electrical   phase   difference   which   is   the 
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negative  of  the  physical  field  crossing  angle  will  result  in  a  gain  of 

V2    at   a  symmetry   point,   regardless   of  coupling. 

Next  the  reciprocal  model  of  SNR  is  compared  with  a 

description   of  reception.      As   indicated   above,   the   general   expression 

for  efficiency  of  an  equivalent  coil  pair  with  coupling     is 
B      [1  +  kAe'',,H"  +  (Ae*  +  ke»)e'*\ 

VR  ^|i  +  kAey,e+p,|2+|(A^  +  ke'V*|2 
Suppose   the   following   situation   applies   for  two   equivalent 

independent    loops: 

S,+N,  ;S,  +  yV,  where   Sa  =  S1e*,(JV1iV2}  =  0  and  ((A'J2)  =  (|W2f). 

Note  that  <|>  here  has  the  inverse  of  its  value  in  transmission  for  the 
same  spin   system.   The  phased   addition  of  the  two   signals  give  the 
SNR 


-  =  V2  ,,  '  .,    ife  =  H».  [4.10] 


If  the  coils  are  coupled  and  then  the  signals   are  phased  and  scaled, 
as  done   for  the   currents  the  result  is 

|S,  +  kS.e'"^  +  A(S.e»  +  kSe*  )e"\ 


|/V,  +kN2e,s> +A(N2  +  k/V,ej1i)e;e| 

S      ll  +  kAe"e+p'  +  (Aeye  +  keM)e»\ 


mj 


[|l  +  kAe/(e+w|2+|Ac"+  ke*f] 


[4.11] 


The      expression   for   gain   obtained   from   reception   and   transmission 

S 

are  exactly  the  same  by  identification  of       B/VRwith      ,    '        .   The 

phasing   and   weighting   which   equalize   current   levels   actually   scale 
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the  signals  according  to  their  levels.     The  effect  of  a  particular 
splitter/combiner    is    different   in   transmission    and    reception    because 
the   coupling   occurs   after  the   power  split  in   transmission   but   before 
the   combination   in   reception.      This   makes   sense   when   one   observes 
that  the  noise  in  each   loop  is  equivalent,   which   means   the  combiner 
adds  SNR  values  which  are  different.     It  is  easy  to  see  that  optimal 
treatment   implies   weighting   them   according   to   their   value.      The 
coefficient  A  does  this  in  the  signal  reception  equation   as  it  equalizes 
current    in    the    transmission    equation. 

So  far  only  the  case  of  equivalent  loops  at  a  symmetry  point 
has  been  considered.     It  can,  however,   be   shown  that  the  conclusion 
that  coupling  does  not  effect  attainable  SNR  is  general  when  there  is 
no    shared    resistance. 

Consider  the  following  situation  for  two  coils  and  fields  at  a 
particular    POI: 

coil  1  has  resistance  R,  and  produces   Bi,   with  input  power 

KfR,- 

coil  2  has  resistance  R2   and  produces  cBi,c'*   with   input 

power    | j,  |"R2. 

B|i.+ci,e*| 
For    this  case  the  efficiency  is  '  "J.  [4.12] 

V|/,|2R,+|/2|:R2 

If  the  power  splitter  is  defined  so  that   i1=Aefii1,      the  proportionality 

Bll  +  cAeVl 
can   rewritten   as    — '  \A  131 

Vr,+a2r, 

Since   the   denominator  has   no  dependence   on   9,  the  expression  can  be 
maximized   by   separately   maximizing   the   numerator   by   choosing   6  = 
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Thus  the  scaled  rotating  field  is  proportional  to 


1  +  cA 

Vr,+a2r,  ' 

cR, 


Through   a     maximization  procedure  one  finds   that      A  =  - — '-  and   the 

R, 


corresponding    efficiency       . —  ^1  +  c2R[/R,  .  When  coupling  is 

A/R, 


^+ke\+c(i2+ke*iJ)e1*\ 

considered,   the  efficiency  is  '  *  .     [4.14] 

y|«,  +  keJpi2|2R,+|i2+  ke'pi,|2R2 
If  i2  =  Aefi\   is  again  defined  and  for  simplicity  of  notation 
K=  ke*  and  a  =  Aefi  are  defined  as  well,     this  can  rewritten  as 


B|]  +  Ka  +  c(a  +  K)(';0| 
V|l  +  Ka|2R, +|a  +  K|2R2 


p 

l  +  c| 

'a  +  K 

]c* 

,  1  +  Ka 

J 

R,+ 

a  +  K 
1  +  Ka 

2 

R2 

[4.15] 


The  form  of  this  expression  is  identical  to  the  expression  without 

a  +  K         R 

coupling  and   by   setting   =  c — -e~*   or     equivalently 

1  +  Ka       R, 


^i    -n 

c — -e 


the  optimum  gain  is     once  again  obtained. 


1-c— Le"J*K 
R, 


Finally,   the  general  case  with   noise  correlation  is  considered,   in 
which  the  fields  at  the  POI  and  resistances  for  two  coils   are  generally 
unequal  and  some  of  the  resistance  is   shared.     The  details  of  the 
derivation   of  noise   correlation   parameters      will   be   investigated   in 
chapter  5.     The  proof  that  coupling  need  not  be  detrimental  follows 
that  above.      First  the  expression   for  rotating   magnetic   field   divided 
by  the  square  root  of  power  is  examined  for  the  case  of  no  coupling. 
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_B_ 


/,  +  ci2e' 


[4.16] 


/|(,|  R,+|i2|  R2+(i, //+!,% )R„ 
The   extrinsic   correlation   which   would   typically   be   associated   with 
the   shared   resistance   is   not   included;   instead   a   hypothetical   ideal 
case   in   which   coupling   has  been   totally  eliminated  is  defined. 
If  i2  =  a/,   is  defined  as  above,  the  following  expression  is 

obtained. 

B  ll  +  cae*! 


VP      yR,+|a|2R2+(a  +  a')Rl2 


[4.17] 


This  can  be  maximized  with  respect  to  a  to  provide  the  best  SNR  (as 
in  Part  II).     When  coupling  is  included  it  is  clear  that 


+  K;2  +  c((,  +Kf,)e' 


B 

Vp 


i,  +  Ki,  |2  R,  +  \i2  +  ki,  |2  R2  +  Rl2  {(/',  +  Ki,  )(i2  +  Ki,  )*  +  (i,  +  k/2  )'(;',+  Ki,  )} 

[4.18] 
If  L  =  ai.   is  defined  it  can  be  observed  that  — ?=  °= 

Vp 

|l  +  Ka  +  c(a  +  K)ey*| 
^|1+  Ka|2R,+|a+  k|2R2+R12{(1+  Ka)(a+  k)*+(1+  Ka)'(a+  k)} 

[4.19] 

and  the  expression  can  be  rewritten  to  show  that 

a  +  K  ^ 
B 

Vp 


+  c 


1  +  KaJ 


R,+ 


O  +  K 


1  +  Kfl 


R2  +  R. 


1  +  Ka     v  1  +  kg 


[4.20] 


.,      .„  ,     a  +  K 

Identification  of  a  and   demonstrates   that   the   form   of  the 

1  +  Ka 

expression   and  therefore   the  maximum  is  the  same  as  the  case  of  no 

coupling. 
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Conclusion 
It  has  been  shown  that  noise  correlation  and  coupling  are 
separable   and   distinct   processes.        Furthermore,   it   has   been 
demonstrated   that     coupling   alone   does   not   represent   any   loss   in   the 
attainable   SNR  if  one  accounts  for  the  parameters  of  the  coupling  and 
uses  the  two  degrees  of  freedom,   weighting   and  phasing,   in   the 
splitter/combiner.      Coupling   is   not  intrinsically   lossy,   rather  it   is   the 
effect  on  the  current  distribution  which  is  of  significance.     Reciprocity 
implies   that     the   current   distribution   gives   all   of  the   information   to 
describe   the   coil.     From  reciprocity  the  desired  end  point  is   known 
and  the  combiner  and  or  the  coils  are  adjusted  to  achieve  this  end. 
The   above   treatment   implies   that   multiple   receiver   systems   can   be 
designed  so  that  all  of  the  effects  of  coupling  are  eliminated   after 
acquisition  providing   the   coupling   constants   are  measured   with   the 
load  in  place.     In  practice  the  smaller  the  coupling  the  easier  the 
compensation,   because  if  the  coupling  is  strong,   a  change  in  the 
characteristics   of  one  coil   affects   all  other  coupled   coils   substantially. 
This   results   in   a  multivariable   solution,   instead  of  a   single   variable 
solution   with   small   perturbation   effects   on   other   coils   which   applies 
for  small  coupling.      In  theory,   however  any  degree  of  coupling  other 
than   unity   can   be   overcome. 


56 


CHAPTER  5 

GENERALIZED  QUADRATURE 

Introduction 

The  first  use  of  an  array  of  coils  to  improve   signal-to-noise  ratio 
(SNR)  in  NMR  appears  to  have  been  by  Chen  et  al.  (32)  with  the  use 
of  a  two  element  radial  array.     This  was  a  so-called  "quadrature" 
volume  coil  in  which  the  individual  saddle  coils  were  oriented  with 
their   fields   perpendicular.      Direct   combination   through   a   90   degree 
combiner  gives  close  to  a     V2   gain  in  SNR  for  this  case.     Several  other 
implementations    of   these   orthogonal    two   element   arrays    have    since 
been  introduced  to  the  literature  (60,61).  The  concept  was  applied  to 
non-volume   geometries   by   Arakawa   et   al.    (33)   with   the    quadrature 
half  saddle  surface  coil.     The  concept  is  the  same  as  for  the  volume 
case  except  that  the  gain  is   approximately    4l   in  a  smaller  region 
where   the   field   lines   are   perpendicular   and   equal   in   magnitude. 
Slight  modification   of  the   geometry   then   led  to  the   completely   flat 
quadrature  surface  coil  (a  two  element  array)  (34).     As  a  result  of  the 
existing   literature,   it  is   commonly   assumed   that   the   benefits   of 
quadrature   reception   exist   only   for   perpendicular   fields   and   with   90 
degree  electrical  phase  angles.     The  validity  of  this  assumption  will 
be  examined   and   the   gain   mechanism   will   be   generalized   to   other 
conditions. 

Recently   arrays   of  several   coils  have   been  used   in   conjunction 
with   parallel   (39)   and   time-multiplexed   (41,   62)   receivers.      In   these 
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applications  the  initial  aim  was  to  extend  the  field  of  view  while 
attempting   to   maintain  the   same   SNR  as   that  obtainable   from  a   single 
element  of  optimal  size  and  position.     The  spine  array  (39)  is  used 
primarily  to  extend  the  field  of  view  (FOV)  in  the  z-direction, 
however,  even   signals  from  coils  far  away  from   a  given  pixel   are 
used  in  forming  the  best  composite  image.     The  volume  abdominal 
array   (63)  consists  of  four  planar  surface  coils,   two  above   and  two 
below  the  body,  which  all,  to  a  great  extent,  "see"  the  same  FOV.     The 
signals  from  all  of  the  elements  are  used  to  obtain  optimum  SNR  gain 
for  a  given  pixel.     Performance  can  be  greatly   increased  with  this 
approach  if  the  coils   are  nearly  independent.     If  two  coils  have  noise 
voltages  which  are  not  exactly  the  same  (within  multiplication  by  a 
complex  constant)   then  the  coils   are  at  least  partially  independent. 
Each  of  the  signals  may  be  separately  acquired  and  after 
reconstruction   the   signals   from   each   pixel   may   be   phased,    weighted, 
and  summed  to  maximize  SNR.     It  can  be  concluded  from  this  work 
that,   in  theory,   the   best  performance  for  a  given  point  of  interest 
(POI)  would  be  obtained  by  receiving  with   as  many   separate  coils   as 
possible  and  combining   their  signals   using   knowledge   of  the  coil 
characteristics.     This  principle  applies  to  fixed  combiners   as  well  as 
independent    receiver    acquisition. 

It    is    important   to   distinguish    between    transverse    and 
longitudinal    arrays    which    are    fundamentally    different   and    for 
optimum   performance   ,   should   be   treated   differently.      Consider,   for 
example,   a  planar  array  of  two  circular  loops  as   shown  below  in 
figure    5-1. 


Figure  5-1.     A  depiction  of  two  circular  surface  coils  assumed  to  be 
independent,    arranged    longitudinally    in   the   direction   of  the    static 
magnetic    field 


In   the   long   wavelength   approximation   the   signals   induced   from 
spins  along  the  bisection  plane,  which  includes  the  z-axis,  in  the  two 
coils  will  always  be  either  in  phase  or  out  of  phase.     The  noise 
correlation   between   the   two   coils   is   independent   of  the   orientation 
with  respect  to  the  main  field.     Optimum  combination  of  the  signals 
will  always  be  with  either  0°  or    180°. 

Consider  a  rotation  of  this  array  by  90°.     When  the  coils  are 
oriented   perpendicular  to  the  main  static   field,   the  phases   between 
signals  originating  from  different  pixels  will  differ  over  a  continuous 
range.     The  noise  correlation  is  the  same  as  in  the  previous  case  and 
optimization  for  each  pixel  in  this  case  will  produce  arbitrary  phases. 
In  general,  this  will  provide  an  increase  in  SNR  over  the  previous 
case.     In  fact,  every  pixel  will  have  a  SNR  which  is  equal  to  or 
superior  to   the  previous  case.     For  the  case  where  optimal 
combination    procedures    aren't    used,    but   instead    the    simple    algorithm 
of  the  square  root  of  the  sum  of  the  squares  of  the  signals  is 
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employed,  there  is  no  difference  due  to  the  above  effects.     There  is 
however    another    difference. 

Consider  the  two  cases  below  in  figure  5-2. 

Z— *>  Z  O 


^_t  "^JS 


Figure  5-2.  A  single  coil  in  the  first  case  (on  the  left)  with  voxel  of 
signal  below  and  displaced  down  the  static  field  direction,  with  the 
spin  orientation  shown.  On  the  right  is  the  same  system  rotated  in 
the   static   field. 


In  the  first  case,  the  magnetic  field  is  never  aligned  with  the 
rotation  of  the  spin  and  the  signal  will  be  much  smaller  than  in  the 
second  case  where  perfect  alignment  with  the  field  will  occur  as   the 
spin  rotates.     This  also  implies  that  the  a  single  coil  has  an  induction 
pattern  which  is   weaker  in  the  z-  direction  than  in  the  x  direction. 
This  stems  from  the  fact  that  the   spins  rotate  in  a  plane 
perpendicular   to   the   static   field.      In   conclusion,   the   transverse   array 
will  be   superior  to  the  longitudinal   array,   assuming   the   coils   in  the 
array   are   the    same. 

If  the  optimum   phasing   and  weighting  of  signals   from   an   array 
is  known  a  priori  for  a  particular  POI,  then   simple  hardware  can  be 
used  to  directly  combine  the  signals.  For  this  POI  the  result  would  in 
theory   be   equivalent   to   a   totally   parallel   acquisition   with   optimum 
phasing  and  weighting.     For  example,  in  the  case  of  an  ideal 
quadrature   volume   coil,   acquisition   of  the   signals   from   the  two 
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perpendicular   coils   through    independent   receivers    followed   by 
optimum  recombination  would  result  in      V2   gain  in  the  center  of  the 
coil  and  would  fall  to  a  gain  of  1  very  near  the  coil  elements.  If  the 
coils    are    combined   with    the    standard   quadrature    combiner   the    result 
is  a/2    gain  at  the  center  but  with  a  more  rapid  fall-off  near  the 
periphery.     This   is  because  both  coils  contribute  equal  noise  for  all 
pixels  but  unequal  signals  in  a  standard  quadrature  combiner.      If  the 
signals    are    independently    acquired    the    appropriate    weighting    can    be 
imposed  to  scale  the  noise  according  to  the  signal  level  for  each  pixel. 

In  chapter  4   it   was   demonstrated  that  coupling   is   not 
necessarily   detrimental,   but      coupling   can   reduce   performance   by 
imbalancing   currents   and  changing   the  phases   of  the   currents.      By 
compensating    for   these   effects    with   the    splitter/combiner   one    can 
prevent  all  loss  for  a  given  POL     Reciprocity  suggests  that  the  current 
distribution    and    the    characteristics    of   the    sample    determine 
everything  there   is  to  know   about  the   system  in  terms  of  SNR   and/or 
magnetic   field.     In  this  chapter  the  effects  of  coupling  on  currents 
will  be  ignored.     Noise  correlation  will  be  defined  as  the  intrinsic 
presence  of  a   shared  resistance   and  the   coupling   which   may   be 
produced   by   this   shared   current   path   will   be   ignored. 

It   seems,   then,   that  all   future   coil   designs   should  consider  these 
findings.      This   work   is   concentrated   on   applying   transverse   array 
technology  to  single  receiver  systems.     Arrays  provide  the  basis  for 
the   best   possible   performance    under   conditions    currently    existing    in 
state  of  the  art  MRI.     The  use  of  arrays  with  fixed  phase  combiner 
networks   will   be   referred   to   as   generalized   quadrature   (GQ) 
combination.    In   the    following   the   importance   of  these   considerations 
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will  be  investigated  and  some  particular  new  coil  geometries  will  be 
discussed.     The  SNR  of  phased  coil  sets  under  various  conditions  will 
be  examined.     First,  the  SNR  of  coil  sets  will  be  calculated  for  the  case 
of  no  noise  correlation.     Since  the  outcome  of  this  is  favorable,  the 
more  common  case  where  the  noise  in  all  coils  in  a  set  is  correlated 
must  also  be  examined.     In  the  last  section  examples  and  new  coil 
designs    are   described. 

Theory 
It  will  be   assumed  in   the  following,   unless  otherwise  noted, 
that   long   wavelength   approximations   apply.      This   implies   that 
radiation   resistance   is   not   part   of  the   central   development   and   that 
the  use  of  the  concept  of  mutual  inductance  is  justified.   It  is  then 
assumed   that  all   noise   comes   from   the   sample,   since   any   other  source 
is,  in  principle,  removable.     The  best  possible  performance  for  a 
given  POI  will  be  obtained  by  the  use  of  as  many  independent  coils 
as  can  reasonably  interact  with  spins  at  that  POI.       To  utilize  this 
approach,     the  degree  of  independence  of  the   coils   and  the  relative 
signals  for  the  coils  induced  from  that  POI  must  be  known.     The 
approach  to  the  calculations  is  straightforward.     The  SNR  will  be 
represented   via   reciprocity    (15)    by   rotating   magnetic   field   per 
square  root  of  input  power  as  in  chapter  4.     In  the  following  the 
attainable   gains   will   be   observed   by  examining   the   magnetic   field 
produced   by   several   coils   when   total   independence   is   assumed.    The 
gain  will  be  defined  as  the  ratio  of  the  efficiency  of     a  combination  of 
similar  coils  to  the  efficiency  of  one  of  the  single  coils.     Next  the  noise 
correlation   issue   will   be   examined   by   calculating   the   change   in   power 
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deposition  when  the   coils   are   placed  on   the   same   loading   sample 

(38). 

Field   Analysis  for  no  Noise  Correlation 

The  case  of  no  intrinsic   noise  correlation  in  multiple  coils   around 
the  same  sample  is  first  considered.     This  is  physically  realistic   in 
many  cases.     For  example,  in  low  field  MRI  and  in  micro-imaging  at 
higher  fields,  the  coil  components  account  for  most  of  the  coil 
resistance.     Alternatively,  one  can  imagine  a  sample  which  is  lossy 
(noise  producing)  only  in  distinct  areas  where  the  total  field  is 
primarily  due  to  a  single  coil  but  the  signal  POI  is  in  a  region  where 
all  coils  can  interact  strongly  with  a  spin.     The  motivation  for  this 
analysis  is  that  gains  greater  than   V2   can  be  obtained  in  theory  for 
certain    physically    realizable    situations. 

If  there  is  no  correlation  of  noise  between  coils,  then  the  coils  do 
not  have   a  shared  resistance.     This   implies  total   independence  of  coils 
which  are  oriented  for  zero  mutual  inductance.     In  each  of  the 
following  three  cases,     power  will  be  supplied  to  the  coil  system  and 
the  resulting  field  which  couples  to  a  spin  will  be  computed.     Recall 
that  magnetic  field  is  a  linear  function  of  source  current.    B  is  the 
rotating  field  at  the  POI  for  unit  current  in  a  particular  coil. 

Consider  first  the  linear  saddle  type  coil  with  resistance  R.     This 
coil  is  depicted  in  Fig.  5-3  from  an  end  view  with  a  vector 
representing   the   magnetic  field  at  the  center  of  the  coil. 
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Figure  5-3.     A  saddle  type  coil  seen  from  an  end  perspective 

The  input  power  produces  a  current  i  in  the  loops  which  produces   at 
the   center  the   field 

B  =  Bie'°"  +  B<e~'™  [5.1] 

where  the  Cartesian  coordinates  have  been  represented  by  real  and 
complex  numbers  and  the  rotating  and  anti-rotating  components 
have  been  explicitly  written.  The  efficiency  for  a  single  coil  is  then 

|Bi1        b 


'lifR      VR 


[5.2] 


Now  consider  the  case  depicted  in  Fig.  5-4  in  which  a  second  saddle- 
type  coil  is  positioned  so  as  to  have  no  mutual   inductance,  resulting 
in  an  angle  <|>  between  the  field  lines  at  the  center  of  the  coil.     As 
expressed  in  chapter  4  it  is  a  matter  of  practical  convenience   to 
minimize  coupling.     The  efficiency  for  this  case  is 
|Bi,+Bi2e*|         B     K+M^l 

VM2R+w2R=^*rw  ■  [5-31 
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Figure  5-4.     Two  reduced  arc  saddle  coils,  oriented  for  zero 
mutual    inductance 


It  should  be  noted  at  this  point  that  the  condition  for  zero  mutual 
inductance   can   be   expressed   as 


f  Br  dS2  =  0. 


[5.4] 


This  is  the     integral  of  the  field  of  coil  one  which  passes  through  the 
loop  defined  by  coil  two.  The  individual  elements  are  assumed  to 
have   the   same   impedance   in   this   configuration   as   they   have 
separately   which   implies   each   of  the   loops   will   have  equal   magnitude 
currents   given   an  equal  power  split.   If  one   of  the  currents  is  delayed 

by  an  angle  6  we  observe  a  gain  of  efficiency  over  the  linear  coil  of 

|l  +  eV'»| 
V2 


[5.5] 
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which  is  clearly  equal  to  V2  if  6  =  -()>. 

Note  that  no  dependence  on  the  value  of  (j>  alone  is  shown.     In 
other  words,  changing  the  physical  angle  of  the  field  crossing  does 
not  change  the  resulting  generalized  quadrature   gain.      For  the  case  of 
no   noise   correlation   quadrature   gain   occurs   even   when   the   fields   are 
not  perpendicular  and  can  occur  even  if  the   field  lines  are  parallel. 
Note  that  <]>   is  variable  under  the   constraint  of  zero  mutual 
inductance,  by  changing  the  arc  length  of  each  coil. 

Now  consider  three  coils  rotated  60  °   from  one  another  as 
shown  in  Fig.  5-5.     For  complete  independence  the  power  will  have 


Figure  5-5.     Three   saddle  type  coils  oriented  for  zero  mutual 
inductance  (at  60  °   rotations)   seen   from   an  end  view   with  vectors 
representing   the  fields   at   the  center  of  the  coils 


increased  by  a  factor  of  3  given  constant  current.     If  the  second  and 
third  coils  have  their  signals  delayed  by  electrical   angles  6  and  i|/ 


respectively,   the  gain  in  efficiency  is 


^3~ 


[5.6] 
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Again  by  delaying  the  currents  so  as  to  exactly  cancel  the 
advancement  due  to  the  relative  position  we  can   force  all  rotating 
components  to  add  in  phase.     Thus  selection  of  9  =  -7i/3  and  \\i  =  -2JC/3 
results  in  a  gain  of  v3. 

The  case  of  4  coils  is  a  straight  forward  extension  and  leads  to  a 
gain  of  2.     This  hypothesis  has  been  tested,  as  will  be  discussed  in  the 
final   section,   and  gains  greater  than   4l   for  multiple  coils  have  been 
found  .  This  is  an  enticing  result  since  in  the  general  case  of  n  coils 
we  see  that  as  long  as  independence  is  assumed  we  obtain   \f  n  gain  in 
SNR  at  the  center.     It  should  be  recalled,  however,  that  the  field 
(which  is  the  basis  for  the  gain)  decreases  in  strength   as   n  increases 
because  of  dwindling  size  of  each  coil,  resulting  in  a  finite  SNR.  It  is 
not  actually  necessary  to  decrease  the  size  of  the  coils,  but  since  the 
current  distribution   is   the   only   important  parameter,   it   is   logical   to 
do  this.     It  should  be  pointed  out  that,  in  this  regime,  gain  can  come 
from   increase   in   the   effective   cross-sectional   area   of  the   conductor 
when  the  resistance  originates  from   ohmic  conductor  losses.      For  two 
coils  the  attainable  gain  is   always  V2    for  equal   strength   fields 
regardless  of  the  physical  angle  of  the  field  crossing.  To  see  if  such 
gains  are  practical,  the  gains  for  the  case  of  correlated  noise  are  next 
calculated. 

Field   Analysis   given  Noise  Correlation 

To   consider  the  case   where   complete  independence   in   the   noise 
is   not  obtained   we  derive   a  relative  expression   for  the   total   power 
deposited  by  the  set  of  coils.     If  there  is  no  correlation,  the  power  is 
simply  the  sum  of  the  individual  powers.     The  cases  are  often 
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specialized  to  that  of  equal  powers  delivered   to  each  coil.   The   degree 
of  correlation  will  be  left  as  a  parameter  in   the   following  discussion, 
since   the   correlation   is   essentially   unpredictable   for   a   typical    sample 
(i.e.   a  human  body)     but  can  be   found     through  direct  measurement 
(39)  (or  measurement  of  the  shared  impedance  (40)).     The   special 
case  discussed  is  the  case  of  equal  correlations  for  all  coil  pairs. 

As   above,   a   certain   time-average   power  is   expended   in   driving 
the  coils.     For  the  identical  coil  sets  above,  the  powers  added  to  form 
the  total  power.       To  derive  the  shared  power  for  a  pair  of  coils,  we 
consider  all   the  power  expended  in  driving  the  coil   as  power 
deposited  in  the  sample.      Several     expressions  for  the  power 
deposited   in   a  sample   by   a   quasi-static   magnetic   field   which   are 
applicable   to   our  problem   have   been   derived   in   the   literature   (38, 
64).     However,  only  a  very  general  expression  is  necessary  to 
accomplish   the  present  goals.     The  average   power  P  expended   in  the 
sample   is   defined   as 

P(B)=  j££,.(BTi)-S/(BTi)dV  [5.7] 

where  the  integral  is  over  the  volume  of  the  sample,  and  the  sum  is 
over   an   arbitrary   number  of  linear  vector   functions   (of  space 
variables)    cjr  BT  is  the  actual  magnetic  field  with  rotating  component 

B.      Linear  functions   have  the  following   useful  properties: 

&1(B1+B1)-$,(B1)  +  $i(B1)  [5.8] 

t(cB,)  =  c^(Bl).  [5.9] 


With   only   these   minimal   requirements   for   the   description   of  the 
power  in   terms  of  the   magnetic   field,   the  relative   values   for  the  total 
power  deposited   by     a  pair  of  coils   which  are   simultaneously   driven 
can  be   obtained. 

P(B,  +B2  )=  |£4,(BIT«,  +B2T;2)-S/(B1TJ,+B2Ti2)<fV  [5.10a] 

=  ll[^Kh)+UKh)}[^{Kh)+UKh)]'dv   [s.iob] 
=  ii%  {Kh  )f + %  {Kh  )f + [§,  (bit/,  )  •  s;(Kh ) + i:(Kh )  •  §,  {Ka  )}<& 

[5.10c] 

=  Jl|^(B1T)f  Kl2  +UK)\2\hf  +[§l(fi„)-§,(fi„)(ft-  +a)]«*v 

[5.10d] 
Since  the  currents   are  not  functions   of  position   they  can   be   removed 
from   the   integral   and   resistances   can   be   defined   from   the   functions 
above,    thus 

P(B,+B2)=  R,|i,|2+R2|^|2+R12  («,//+«,%)  [5.11] 

where    R,  =  JX|^(b,t)|W,R2  =  JS^B^lW  and 

Rl2  =  JX^:(Bit)'^,(B2T)^V.  It     should  be  clear  that  the  resistances   R, 

and  R2  are  the  isolated  coil  resistances  in  the  presence  of  the  sample 
but  with  the  other  coil  absent.  Note  that  the  condition  of  zero  mutual 
inductance   (eq.    [5.4]   )  does   not   necessarily  imply   that  this   shared 
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power  term   is   zero.     This   derivation   also  demonstrates   that   the   field 
functions  which  produce  the   shared  power  can  he  orthogonal  in 
space  and/or  in  time.      By   electrically  phasing  them  different  by  90° 
orthogonality   is   obtained   regardless   of  the   spatial   overlap   of  the 
fields   .   If  two  adjacent  loops  are  combined   so  that  their  currents  are 
in  phase,   it  is  apparent  that  the  combination  looks  much  like  a  loop 
twice   as   large  with   a  corresponding  increase   in  depth   penetration 
over  a  single  coil.     If  the  loops  are  driven  with  opposite  phases  then 
the   net  field   turns   back  on  itself  resulting   in   less   depth   penetration 
than  a  single  coil.  The  power  expression  above  makes  these  effects 
explicit  and  describes  the  net  effect  of  the  vector  addition  of  the  field 
in  terms  of  noise  equivalent  resistance.     From  the  form  of  the  power 
expressions  it  is  clear  that  by  the  triangle  inequality,  Ri  +  Ri  >  2Ri2- 
This   ensures   that   the   total   noise   resistance   cannot   be   negative. 

Although   the   resistance   discussed   here   is   from   the   inductive 
loading  from   the   sample  it  should  be  pointed  out  that  all  effective 
losses  will   behave  the  same  since  all   are  representable  as  quadratic 
forms  of  current.     Ohmic   losses  in  separate  probes  have   no   intrinsic 
correlation   however   orthogonal   modes    which    occupy   the    same 
physical  conductors  in  general  will.     Typical  implementations  of 
these,   such   as  the  birdcage,   phase  orthogonal   modes  by  90°    which 
eliminates    the   effect    from    observation. 

The  calculation  of  the  efficiency  for  two  coils  (eq.  [5.3])  is 
repeated,  this  time  including  the  noise  correlation.      In   the  following  it 
is   assumed   that    i2  =  \e*   and  isolated  resistances  are  equal  (R]  = 
R2=R).  This  implies  that 
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B__ B|/,+/,cVa| 

/P      ^R|i,|2+R|/>i8|:+R1,(|j|fc-JB+|/1|2cje) 

Bll  +  eVl 


/Ri/2(l  +  acos6) 


[5.12] 


o 

where    a  =  — —  and  will  be  called  the  correlation  coefficient.     In 
R 

general  it  is  not  true  that  optimum  SNR  is  obtained  by   setting  the 

electrical  angle  equal  to  the  negative  of  the  field  crossing  angle,   since 

now  a  more  complex  function  of  electrical  angle  6  is  obtained.     From 

eq  5.12,  it  appears  that  in  all  two  element  cases  the  possible  gain 

approaches   infinity   as   the  correlation   coefficient   a  approaches     1, 

unless  <(>  is  zero.  Of  course,  for  a  truly  to  reach   1   the  fields  would  have 

to  be  exactly  the  same  and  if  is  necessarily  zero  for  all  locations. 

Therefore  while  large  gains  may  be  realized  there  is  a  real  limit  on 

the   gain. 

For  the  general  case  of  three  coils  the  power  calculation  done  for 

two  coils  can  be  easily  modified  with  the  result  that  the  total  power 

is 

|zf  [Rl  +R2  +  R3+ 2cosGRi2  +   2cos(y-6)R23   +  2cosv|/R]3].     [5.13] 

If  we  again  make  the  assumption  that  Ri  =  R2  =  R3  =  R  and  further 

assume   that   the   magnitudes   of  the   shared   resistances   are   equal   (R)  2 

=  R13  =  R23)  ,  as  is  required  for  the  worst  case  analysis  of  total 

correlation,  we  find  again  that  a   =   1   for  total  noise  correlation  where 

a  is  defined  as  R12/R.     Note  that  the  total  power  can  mathematically 

reach   the   value   9|('|2R  when  all  electrical  angles  are  chosen  as  0°  and 


for  total  correlation  (a  =  1).     This  case  implies  a  maximum  gain  of  one 
for   the   system. 

For  the  three  element  volume  coil  in  figure  5-5,  it  is  not 
possible  to  obtain   this  mathematical  maximum  for  the  noise.      For  this 
case,  the  defined  orientation  for  correlation  from  coil    1   to  coil   3  is 
opposite  that  of  1   to  2  and  2  to  3,  thus  we  obtain  the  total  power 

|(f  R[3  +  2<x(cos9     +   cosO-6)   -   cosy)].  [5.14] 

For  the  same  phase  angles  chosen  in  the  case  of  no  correlation  (9  =  - 
7T./3,  y;  =  -  2TC/3)  we  obtain  a  net  gain  of  -^3/2    for  the  worst  case  of 
total    correlation   (a  =  1).     It  is  interesting  to  note  that  even  for  total 
correlation  one  obtains  a  GQ  gain  relative  to  a  single  coil  element.     For 
arbitrary  correlation,   the   total   power  is      |/|23R[1  +  a]  for  the  above 
choice  of  phase  angles  and  is  the  largest  attainable  power  for  a  given 
current  magnitude,  but  as  we  will  show  in  a  later  section  this  choice 
produces  the  highest  SNR  in  the  center  as  well. 

There   has   been   considerable   debate   about   noise   correlation   in 
relation   to  mutual   inductance   (37,38).        Coupling   and  correlation  have 
been  treated  as  distinct  and   not  necessarily  related.   Some  empirical 
evidence     suggests  that  it  is  impossible  to  have  zero  crosstalk  and 
still   have   noise   correlation.      It  has   been   recently   claimed   that  the 
presence   of  noise   correlation   between   coils   necessarily   implies   that  a 
signal  path  exists  between  these  coils   (40).     This  was  disproved  in 
chapter  4.     Furthermore  it  was   shown  in  chapter  4   that  the  effects  of 
coupling  can  be  eliminated  in  other  ways.     Maximal   isolations   for 
overlapped  planar  coils  on  the   order  of   10   to   20  dB   have  frequently 
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been    observed   for   sample    dominant   loading    conditions    where    both 
coils  are  equally  loaded  by  a  single  large  sample.     For  sample 
dominant   loading   in   which   two   separate   samples   were   used   to   load 
the   coils   independently,   the   attainable   isolation   was   much   greater. 
These  results  were  obtained  using  two  loops,  each  on  a  different 
substrate.     The   isolation   between   coils   is   measured  as   the  degree  of 
overlap   is   varied.        Furthermore  the   author  has  found  that  the  zero 
mutual   inductance   overlap   for   the   unloaded   system   does    not   give 
maximum   isolation  under  loaded  conditions.     It  is  clear,   then,   that  a 
significant  amount  of  signal  may  be   transferred  from  one   coil   to 
another  when  we  have  the  condition  of  substantial   noise   correlation. 
In  the  next  section  we  show  real  examples  in  which  excellent  gains 
are   obtained   in   the   presence   of  substantial   coupling   caused   by   the 
shared    resistance    of   noise    correlation    supporting    the   theoretical 
developments  in  chapters  4   and  5. 

Applications    and    Examples 

Three   Phase  Volume  Coil 

It  is  of  interest  to  know  how  the  new  three  phase  design 
compares   in   performance   to   standard   designs.      The   techniques 
derived  above  will  be  used  to  obtain  approximate  limits  for 
performance  comparisons.     Coupling  will  be  ignored  for  this  example. 
It  is  assumed  that  the  sample  accounts  for  essentially  all  of  the  noise 
but  that  the  sample  is   not  extremely  close   to  the  coil  elements.  Under 
these   conditions   we   conclude   that   a   standard    120°  arc   saddle   coil   can 
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be  approximated  by     two  60°  arc  saddle  type  coils  which  are  slightly 
overlapped   for  minimal   coupling.      This   produces   approximately   a   60° 
physical   angle   (<|>)  and  an  electrical  phase  difference  (6)  equal  to  0° 
completes   the   construction.      Since   the   currents   are   in   phase   the 
elements   which   overlap   will   essentially   cancel   each   other   in   the 
region   of  interest   and   the   result   is   approximately   an   equivalent 
larger  saddle  coil.     The  three  phase  coil  (Fig.  5-5)  is  shown  in 
comparison   with   the   same   electrical   phase   angles   discussed   above. 
From  eq.  5.5  we  see  that  the  gain  is  a/3/2   compared  to  the  60°  arc 

saddle.     The   3-phase  coil  (figure  5-5)   with  optimal  phase   angles 
produces    a/3  gain  (see  eq.  5.6)  relative  to  a  single  60°  arc   saddle   coil. 
Multiplication   of  the    120°  arc  linear  saddle  case  by   a/2    shows  we 
obtain   exactly  the   same   performance   in   the   center  for   the   standard 
quadrature  saddle  coil  and  the   3   phase  coil.     These  results  are 
summarized  in  Table  5-1.     Although  the  analysis  is  not  exact  it  is 
clear  that  the  greater  gain  of  the  GQ  coil  may  not  correspond  to 
improved   SNR   over  other  designs. 

Four  Phase  Volume  Coil 

Consider  the  case  where  two  saddle  type  coils  with  reduced  arc 
lengths   are  oriented  for  minimal  coupling   at  a  45°   physical  angle  for 
<t>.     This  coil  is  depicted  in  Fig.  5-4.     Again  assume  that  a  loading 
sample   accounts  for  essentially  all   of  the   noise  but  that  there  is 
substantial    spacing   between   the   coil   elements   and   the   sample. 
Further,   take   Ri    =    R2.       An  estimate  is  obtained  for  the  increase  in 
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SNR  of  the  optimally  phased  pair  and  a  single  saddle-type  coil  with 
arcs  which  are  approximately  twice  as  long  or  about  90°. 

Table  5-1.     Approximate  GQ  gain  comparison  of  a  standard  linear 
saddle  coil,   a  quadrature  saddle  coil,   and  a  three  phase   volume  coil 

all  with  a  60°   arc   saddle  type  coil  as  reference. 

independent  coil  case  correlated  noise  case 

•    ,    .  /TI 7w*      fin  Vl+cos(60°)       ^3/2 

equivalent  Vl  +  cos(60  )  =  a/3/2  - ,         =    , 

linear  saddle     (120  °) 


equivalent 

quadrature   saddle   (120   °)     V3 

3  phase  (60  °)  V3 


V3 
Vl  +  cc 

V3 


vr+ 


(/. 


Table  5-2  indicates  the  GQ  gains  over  a  single  45°  arc  saddle  type 
coil. 

Consider  the  phased  pair  of  45°   arc  saddle  type  coils  for  a 
typical  correlation  coefficient.     For     a  =  0.3,  the  GQ  gain  maximum 
occurs  at  about  9  =  -  75°  and  is  about  1.32.     If  it  is  assumed  that  the 
sample  prohibits  coverage  of  more  of  the  azimuthal  angle  than  a  90° 
saddle  coil,  then  a  gain  of  15%     in  SNR  is  obtained  by  using  the 
phased   pair  since   the   linear   saddle   with   0.3   correlation   coefficient 
produces  a     relative  gain  of  1.15. 
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Table  5-2.     Approximate  GQ  gain  comparison  of  a  linear  90°  arc 
saddle  type  coil,  a  phased  pair  of  45°   arc   saddle  coils,   a  quadrature 
90°  arc  saddle  type  coil,  and  a  four  phase  coil  all  referenced  to  a 
single   45°  arc  saddle  type  coil. 

independent  coil  case  noise  correlation  case 

equivalent 

linear   saddle    (90°)  Vl  +  cos(45°)  =  1.307  V   +cos(_ — )_=       

VI  +  a  VI  +  a 


phased  pair  V2  (8  =  -  45°) 


equivalent     quadrature 


(90°)  saddle  J  2  Vl  +  cos(45°)  V"2 


4  phase  coil 


Vl  +  cos(45°+9) 
Vl  +  acos(0) 

Vl  +  cos(45°)  _     1.85 
Vl  +  a  Vl  +  a 


Vl  +  acos(45°) 


A  perpendicular  pair  can  be  added  to  these  coils  to  obtain  another 
■J  2  increase  in  SNR  by  phasing  at  90°.     This  requires  selection  of  8  =  - 
45°  for  the  phased  pair  to  provide  uniformity.     The  GQ  gains  are 
indicated  in  Table  5-2.     From  this  analysis  it  is  concluded  that  the  4 
phase  coil  has  improved  SNR  compared   to   a  quadrature   (90°)   saddle 
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coil  in  the  range  of  8%  (a  =  0)  to  17%  (a  =  1).    Furthermore  the 
homogeneity  and  SNR  of  the  4  phase  coil  is  essentially  that  of  a 
quadrature    8   element   birdcage    (15)    since    the   current   distribution   is 
approximately    the    same. 

These  principles   are  tested  by  using  a  4  element  array  each   of 
which  was  a  curved  surface  coil  with  about  a  50°  arc.     This  is  one 
half  of  the  cylinder  of  coils  depicted  in  fig.  5-6. 


Figure  5-6.     An  eight  coil  array  seen  from  the  end,  with  each  element 
rotated    45°    with   respect  to   the   previous   element 


The  Q's  of  these  coils  were  reduced  to  about  4  (from  about  60)  by 
adding  resistors  in  parallel  with  capacitors.     This  made  the  isolation 
20  dB  or  better  for  each  pair.     The  coils  were  placed  around  a 
nonloading   sample   and   an  image   was   made   using   a  4-way   combiner 
which   added  their  signals   at  45°  increments.      Images   were   also 
produced  with  each  individual  coil  with  the  other  coils  in  place  and 
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still  connected  through  the  combiner.     The  two  empty  ports  of  the 
combiner  network  were  terminated  with  50  £2  loads.   Because  of  the 
resistors  this  is  a  case  of  essentially  no  correlation  of  noise.     Fig.  5-7 
shows  each  of  the  individual  coils  images,  the  combination  of  all  coils 
in  rotating   mode  and   in  anti-rotating  mode.      Table   5-3   shows 
isolations   and  table   5-4   shows   SNR  values. 

Table  5-3.      Isolation  measurements   (dB)  of  4  element  half-cylindrical 
array  under  loading  condition  used  in  images  of  fig.   5-7. 


Coil  1 

Coil  2 

Coil  3 

Coil  4 

Coil  1 

- 

20 

20 

20 

Coil  2 

2  0 

- 

23 

20 

Coil  3 

2  0 

23 

- 

21 

Coil  4 

2  0 

20 

21 

_ 

Table  5-4.         SNR  values  for  4  phase  coil  measure  at  position  L41, 
A70  for  all  images  of  fig.  5-7 

SNR  Gain   of  rotating   mode   over   individual   element 

Coil  1  90.9  1.95 

Coil  2  110.1  1.61 

Coil  3  95.5  1.85 

Coil  4  109.7  1.61 

This  same  concept  was  applied  to  an   8-phase  coil  with 
elements  as  in  fig.  5-6.     In  this  case,  however,  local  samples  (7  M 
sodium   chloride)   were   produced  for  each   coil   to   preferentially   load 


Figure  5-7.     Top  four  images  from  individual  coils- 
bottom   left   GQ   addition,   bottom   right   anti-phase 
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only   it  and  reduce  correlation  of  noise  between  coils.      Isolation   was 
not  as  high  as  for  the  4  element  case  nor  was  reduction  of  correlation 
as  nearly  complete.     Fig.  5-8  shows  the  8-way  GQ  addition  of  all  the 
loops  in  fig.   5-6,   while  figures   5-9a  and   5-9b   show  the   images  made 
by  each  individual   coil.     Tables   5-5   and  5-6   show   relevant  data. 

Table  5-5.        Isolation  measurements  for  8  element  array  loaded  by 
sample  used  in  images  of  fig.  5-8 

Coil  1  Coil  2      Coil  3      Coil  4      Coil  5      Coil  6      Coil  7     Coil  8 

Coil  1        -  26.0 

Coil  2  26 

Coil  3  18.3  22.6 

Coil  4  22.4  18.5 

Coil  5  25  23.0 

Coil  6  21.7  25.2 

Coil  7  17.7  23.2 

Coil  8  24.7  17.8 

It  has  been  shown  with  these  examples  that  the  assertion  of  multiple 
coils  producing   greater  than    V2   gain  in  the  case  of  no  noise 
correlation   is   valid  and  in  the  next  section   the  claims   for  correlated 
noise    are    addressed. 


18.3 

22.4 

25.0 

21.7 

17.7 

24.7 

22.6 

18.5 

23.0 

25.2 

23.2 

17.8 

- 

23.4 

17.7 

27.8 

25.2 

22.3 

23.4 

- 

28.0 

16.9 

21.8 

23.8 

17.7 

28.0 

- 

25.0 

17.9 

22.2 

27.8 

16.9 

25.0 

- 

26.7 

17.1 

25.2 

27.8 

17.9 

26.7 

- 

22.0 

22.3 

23.8 

22.2 

17.1 

22.0 

. 

so 


Figure  5-8.  GQ. combination  of  eight  coils 


Figure  5-9a.  Individual  images  from  four  coils 
used  with  coils  of  figure  5-9b  to  make 
figure  5-8 


S2 


Figure  5-9b.  Individual  images  from  four  coils 
used  with  coils  of  figure  5-9a  to  make 
figure  5-8 


Table  5-6.     SNR  measurements  made  at  location  R38.  P38  in  images 
of  fig.   5-8.     Corrected  values  use  extra  losses  observed  in   network 
utilized   for   combination. 


raw  SNR 

Coil  1 

255 

Coil  2 

214 

Coil  3 

248 

Coil  4 

258 

Coil  5 

255 

Coil  6 

260 

Coil  7 

267 

Coil  8 

207 

8-phase 

420 

Surface 

Coils 

Raw  Gain        Corrected  SNR      Corrected  Gain 


1.65 
1.96 
1.70 
1.63 
1.65 
1.62 
1.57 
2.03 


255 
214 
228 
212 
246 
231 
225 
185 
420 


1.65 
1.96 
1.84 
1.98 
1.71 
1.82 
1.87 
2.27 


All  of  the  examples  of  GQ  coils  given  so  far  have  been  volume 
coils,  but  the  two  element  coil  case  is  applicable  and  most  useful  for 
surface  coils.     Consider  Fig.   5-10,  which  shows  the  directions  of  the 
B  l    fields   (at  a  certain  location)  produced  by   two  overlapping   surface 
coils.     It  is  obvious  that  the  physical  angle  (])  changes  as  the  POI 
moves  along  the  symmetry  axis.     This  is  confirmed  by  fig.   5-11 
which   shows   a  sequence   of  images   obtained  by  changing   the 
combination   phase      to   various    antiquadrature   values.    Note   that   the 
location  of  cancellation  moves  along  the  central   axis. 


Figure   5-10.      Representative   field   lines   from   two   overlapped   surface 
coils  crossing  at  angle  <)>  on  the  symmetry  axis 


Equation  5-12  implies  that  for  regions  where  the  angle  0    approaches 
180°   the  gain  can  be  greater  than   a/2    and  where  the  fields  cross   at 
angles  that  approach     0°   the  maximum  gain  will  be  less  than   a/2. 
This  further  implies  that  one  can  choose  the  phase  angle  which 
produces  the  optimum  SNR  for  a  given  depth.     This  conclusion  was 
experimentally   verified   with   a  pair  of  identical   surface   coils   with 
strong  noise  correlation.     By  phasing  the  pair  at   145°  we  were  able  to 
obtain  a  raw  gain  in  SNR  of  59%.     Fig.     5-12  shows  these  three 
images,      table  5-7      shows  relevant  parameters  for  these  coils,   and 
Table   5-8   shows   actual   and   predicted  gains   for  various   cases. 
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Figure   5-11.      A   sequence  of  images  produced   with   different 
combiner  angles   indicating  the  different  crossing   angles   at  different 
locations 


8  6 


|i|P^£^||:*|iJ::.'::':'''':: 

XtteZ 

T:;  ., 

■ 

:  M^i? 

* 

MH: 


— 


V:! 


Figure  5-12.     Surface  coil  images  —  top:  right  coil;  middle:  left  coil; 
bottom:      145°   combination  of  both  coils 


Table  5-7.     Bench  measurements  of  individual  coils  and  phased  pair 
under  conditions  of  magnet  test  from  images  of  fig.  5-11. 


individual      coils: 

Resistance  of  untuned  loop  measured  with  vector  impedance 

meter 

Rr  =  R2  =  18  n 
Relative  current  measured  in  loop  with  network  analyzer  and 
magnetic    field    probe 

Ij  =  I2  =  -  28.4  dB 


pair    of    coils: 

isolation  before  combination  =  12.5  dB  =>    k  =  0.237,    (3  =  180° 
Resistance  of  series  aiding  (0  °)  combination    R=  49  £2 
Resistance  of  series  opposing  (180  °)  combination     R=  21  Q 

Currents  measured  in  loops  with  145°   phase   difference   with 
same  power  as  for   individual  case 

II  =  -28.7  dB      12  =  -  29.1  dB 


Table  5-8a.     Actual  SNR  measured  at  R41   P58  for  coils  of  fig.  5-12. 

SNR 

Coil  1  21.5 

Coil  2  21.3 

Combined  coil    145°         34.1 


Table   5-8b   Actual  versus  predicted  gains 

Actual  Gain  159 

Maximum  gain  predicted  from  current  measurements  1.89 

Maximum  gain  predicted  from  theoretical  considerations  1.71 

In   this   case  three  different  coils   were   constructed  on  a 
standard  spine  coil  form  from  MRI  Devices  Corporation  of  Waukesha, 
WI.     First  only  the  left  loop  was  tuned,  in  the  next  case  another  board 
had  only  the  right  loop  tuned  and  finally  on  another  board  both  loops 
were  tuned  and  combined  so  that  their  currents  differed  by      145°. 
The   tests   were  performed   with   an  open  positioning  table,   so  that  the 
placement   of  the   coils   and   sample   were   extremely   reproducible.   The 
loader  was  a  Plexiglas  box  of  saline  and  copper  sulfate  of  dimensions 
6.5  cm  x  35.5  cm  x  46  cm  and  was  positioned  so  that  the  saline  was 
about  3  cm  from  and  centered  above  the  coils.  The  condition  that  Ri 
=  R2  =  R  was  therefore  satisfied,  so  that  eq  5.12  was  used  to  compute 
the  expected  gains  from  this   system.   The  correlation  coefficient  a 
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was  computed  from  the  value  of  resistance  measured,  giving  a  = 
0.39.  The  actual  gain  was  lower  than  the  predicted  gain  because  of 
small  combiner  losses  and  because  the  field  crossing  conditions  for 
maximum  gain  to  be  realized  do  not  precisely  occur  in  our  coil 
design.   The   agreement   of   the    current   measurement   (which    includes 
combiner  losses)  gain  and  theoretical  gain  is  fair. 

Summary   of   Attainable   Gains 

Below  is  a  description  of  various  domains  which  differ  in  terms 
of  attainable  SNR.     Descriptions  are  based  on  a  transmission  model. 

1.  No  power  deposited  in  sample.  Ohmic  coil  losses  only. 

The  total  power  expended  in  driving  the  coils  is  the  sum  of  the 
individual  ohmic  power  losses  in  each  coil.     Reciprocally  there  is  no 
correlation  of  noise.     Therefore  if  n  equivalent  coils  can  be  uniformly 
positioned   around   a   sample   with   correctable   coupling   one   can   obtain 
Vn   gain  in  SNR  over  the  unit  coil. 

2.  Some  power  deposited  in  sample. 

The  amount  of  gain  depends  on  specific  conditions.     Noise 
correlations  will   never  be  total   because   some  power  goes   into   ohmic 
losses.   Some   gain   is   acheivable  for  any  relative  orientation. 

3.  All  power  deposited  in  sample. 
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a.  non-uniform  fields  in  sample 

For  non-uniform  fields  the  gains  again  depend  on  specific  conditions. 
The  degree  of  correlation  is  arbitrary  and  since  the  fields  are  non- 
uniform it  is  possible  to  obtain  large  gains  if  the  dominant  fields  in 
the  sample  are  much  different  in  crossing  angle  than  the  angle  at  the 
POL 

b.  uniform  fields  in  sample 

This  narrow  domain  is  the  common  ideal  case  for  volume  coils  and  is 

examined   in   detail. 

Uniformity  of  the  magnetic  fields  allows  removal  of  the   BT 

field  terms  from  the   spatial  volume  integral.     This  results  in 

modification  of  eq  5.10d  to 


P(B, +B2): 


|bit|2|/,|2  +|b2T|2  Up  +b1t  -b2T(m' +ft)]jlfc(i)|>  • 

[5.15] 


For  BIT  and  B2T  equal  except  for  a  rotational  change  in  position  by 

angle  <j>  we  obtain  a    correlation  coefficient  of  a  =  cosifi.  Equation  5.12 

becomes 

B  Bll  +  eVM 


VP        VR^2(l  +  cos(j)cosej' 


[5.16] 
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This  expression  is  maximized  by     9  =  -ty  and  the  gain  is      —       —  ■ 

For  0=0  the  coils  noises  are  totally  correlated  which  results  in  no  gain, 
whereas  for  <|>=7i/2,  the  attainable  gain  is   a/2.     Furthermore  for  the  n- 
phase  volume  coil  in  which  we  continue  the  assumption  of  all  power 
to  the  sample  and  uniform  fields  throughout  the  load,  we  find  that 
the  maximum  gain  remains  fixed  at   V2  for  n>2  and  cannot  be 
increased. 

Conclusion 

It  was  demonstrated  that  standard  quadrature  reception  is  a 
special  case  of  a  more  general  principle  which  describes  weighting 
and  phasing  for  several  coils  to  obtain  optimum  SNR  at  a  given 
location.     This  principle  was  demonstrated  for     POI's  which  were  all 
at  locations   where  the  coils  all  provided  approximately  equal 
magnitudes  of  signal.  The  method  could  easily  be  extended  to  other 
regions.     The  concept  of  independence  of  acquisitions  led  to  the 
consideration  of  choosing   a  particular  phase  and  weight  for  each  coil 
in  an  array  and  then  combining  the  signals  into  a  single  channel  data 
stream.     The  result  is  optimum  performance  at  some  specific  location 
and   less   than   optimum   performance   elsewhere.      The   same   array 
used   in   parallel   acquisition  provides   the   ability   to   perform   the 
combining  in  software  on  a  pixel  by  pixel  basis  leading  to  optimal 
performance   everywhere.    The    number   of   channels    is   limited 
however,   and  it  may  be  that  each  element  of  a  parallel  acquisition 
should  be  a  GQ  array. 
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The    approximate   theoretical   performance   of  three    and   four 
phase  cylindrical  volume  coils  was  examined  and  it  was  found  that 
the   four  phase  design   appears  to  provide  both  excellent  homogeneity 
and  SNR.     Theoretically,  the  cylindrical  geometry  allows  for  canceling 
of  mutual  inductance;  however,  it  was  previously  shown  that 
excellent   performance   can   be   obtained   even   with   significant   coupling 
if  compensation  is  made  for  the  phase  and  magnitude  of  the  coupled 
signal.     This  implies  that  arrays  such  as  the  abdominal  array 
described  by  Hayes  et  al.  (63)  can  be  quite  effective  even  though  the 
elements  appear  to  be  coupled.     Isolations  on  the  order  of  10  to  20 
dB  have  proved  to  be  sufficient  for  good  performance  on  real  GQ 
coils.     The  correlation  coefficient  was  derived  by  the  use  of  long 
wavelength   approximations   that  may   not  be  justified   at  frequencies 
substantially  higher  than  64  MHz  which  is  currently  the  highest 
frequency  commonly  used  in  clinical  MRI.  The  usage  of  mutual 
inductance   also   must   be   assessed   at  higher   frequencies.   This   method 
is  also  helpful  for  comparisons  of  standard  designs. 

The  two  element  coil   analysis   (eq.   5.12)   demonstrated   that  a 
pole  in  SNR  occurs  in  the  case  of  total  noise  correlation.    This  leads  to 
potentially  high  SNR  gain  even  for  coils  which  have  very  strong  noise 
correlation  if  the  field  lines  at  the  POI  do  not  differ  by  exactly  the 
phase   of  correlation   (0°  or   180°).     Furthermore,  for  a  pair  of  surface 
coils,   a  depth  dependent  variation  of  geometric   phase   angle  can  be 
produced.     This  suggests  that  a     variable  phasing  system  could  be 
produced   which   would   allow   optimization   at   a  particular  depth.      It 
might  also  be  useful,  in   some  instances,   to  consider  transmission   with 
one   phase   and  reception   with   an   unrelated   phase   to   produce   a  given 
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excitation  profile  and  possibly  a  different  reception  profile.  This 
could  prove  useful  for  some  kinds  of  surface  coil  localization  and 
subcutaneous    fat    suppression. 


CHAPTER  6 

INDEPENDENT  SIGNAL  ACQUISITION 

Introduction 

The  previous  3  chapters  deal  with  optimal  combination  of 
signals  into  a  single  channel.     It  was  known  that  each  coil  received 
signal  from  various  regions  of  the  sample  and  noise  from  similar 
regions.  The  relative  phase  and  amplitude  of  the  signal  acquired  by 
each   coil  was   used   with  overall   correlation   parameters   to   define   the 
single   combination   which   produced   the   highest  SNR  at   a   particular 
location.     This  choice  of  a  single  point  for  optimization  is  restrictive; 
however,   there  historically  has  always  been  a  trade-off  between  SNR 
and  field  of  view.     Note  that  the  philosophy  of  knowledge  of  the 
problem   was  highly  relevant  to  the   solution.      Optimal   combination 
required   an   approximate   knowledge   of  the   coil   field   characteristics, 
crosstalk   and   intrinsic   noise   correlation. 

This  chapter  deals  with  increasing  the  knowledge   available  for 
reconstruction  of  the  image  by  tagging  the  signals  from  each  probe  in 
a  recognizable  way,  so  that  they  can,  in  software  reconstruction,  be 
optimally  added  for  all  locations  in  the  resultant  image.     This  is 
theoretically   equivalent   to   an   optimum    combiner   designed   for   each 
POL     Roemer  et  al.  (39)  were  pioneers  in  this  area  and  their  work 
resulted  in  the   most  general   method  for  this   process.     Their  system 
employs   a  separate  receiver  for  each  coil.   Obviously   the   signals   from 
each   coil   remain  distinct,   allowing   a   separate  reconstruction   for  each 
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coil  and  then  a  pixel  by  pixel  optimized  addition  of  the  images  based 
on  the  SNR  in  each  pixel  of  each  image.     There  are,  however,  many 
ways  to  label  each  signal  so  that  the  signals  are  effectively 
orthogonal  and  can  be  combined  in  a  single  channel  receiver  and 
then  be  decomposed  into  their  respective  signals  and  finally 
recombined  to  produce  the  final  image.     Several  methods  will  be 
discussed   in   this   chapter. 

Time    division    multiplexing 

The  method  of  Roemer  et  al.   (39)  of  using  multiple  receivers  to 
obtain  information  from  multiple  coils  without  ambiguity  is  a  classic 
engineering   example   of  a   parallel   implementation.      The   standard 
alternative  to   this  is  the  serial   implementation  which   will   be   called 
time  division  multiplexing  (TDM).     The  idea  is  to     break  up  the 
signals  into  different  pieces  in  the  time  domain  thus  make  an 
alternating  series  out  of  the  signals'  time  pieces.     Since  the  signals  are 
later  sampled  it  is  only  necessary  to  maintain  each  signals  integrity 
at  the  sampling  time.     If  the  sample  and  hold  circuitry  is  fast  enough, 
multiple  signals  can  be  forced  into  the  same  total  time  as  is 
necessary  for  one  of  the   signals.      Other  requirements  will  be  derived 
below. 

Switching  in  TDM  can  be  modeled  as  the  multiplication  of  a 
square  wave  of  height   1   and  duty  cycle   1/N  (for  N  coils)  with  the 
signal  of  interest.  This  is  done  at  rf;  however,  all  of  the  mixing  to 
lower   frequencies    and    quadrature    detection    all    involve 
multiplications  as  well  (by   sinusoids).     Since  multiplication  is  a 
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commutative  process,  the  switching  can  be  analyzed  as  if  it  occurred 
at  the  lowest  frequency  of  the  receiver  system.  This  applies  as  long 
as  there  are  no  filters  which  have  bandwidths  on  the  order  of  the 
switching  rate.     In  general  this  is  true  (bandwidths  about   1-2  MHz) 
until  the  last  audio  filter.     Figure  6-1   demonstrates  a  block  diagram 
of  a  multiplexed   receiver. 


Switch 


standard 


receiver 


Coils    Preamplifiers   Bandpass 

filters 

Figure   6-1.      A  block  diagram  of  a  time-multiplexed  receiver 


The  simplest  case  of  TDM  is  the  switching  of  two  bandlimited 
audio  signals  as  shown  in  figure  6-2.  In  frequency  space,  the 
multiplication  by  a  square  wave  spreads  the  information  out  over  all 
frequencies.  This  is  depicted  in  figure  6-3  by  demonstrating  the 
magnitude  of  the  Fourier  transforms  of  the  signals  of  interest.  For 
audio,   f0  equals  zero.     The  process  of  the  spreading  of  information  to 
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Figure    6-2.        Mathematical    equivalent    of    switching    between    two 
signals   and   a   schematic   representation   of  this   process. 


all  frequencies  is  indicative  of  the  existence  of  a  related  inverse 
phenomenon. 

Any  signals  outside  the  band  of  interest  may  be  folded  into  the 
frequency  band  of  interest.     Since  the  signals   are  multiplied  by  a 
switching  function  at  rf,  it  is  also  necessary  to  band  limit  at  rf,  so  that 
noise  from  outside  the  band  will  not  be  folded  in.     There  is  a 
requirement   then   of  extremely   narrow   band   filters    which   have   the 
same   bandwidth   as  the  rf  signal   and  have   steep  fall-off  skirts. 
Fortunately   crystal  filters  have  these  properties  as   is   obtained  from 
the   transmission   profile  of  a   standard  crystal  filter  as  discussed  in  a 
later  section.     A  3  to  4  dB  loss  in  the  pass  band  is  negligible  when  the 
signals   are  first  subjected  to  preamplification.     This  gain  is   necessary 
so  that  the  noise  is  treated  like  a  signal  by  the  filter  and  therefore 
attenuated. 
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Figure   6-3.      Magnitudes   of  Fourier   transforms   of  the   original    signal 
(top)   and   the   switched   signal   (bottom) 


If  the   signals  depicted  above  were  preamplified  and  filtered  at 
RF  then  the  SNR  of  each  image  in  frequency  space  is  almost  the  initial 
SNR.     If  the  information  in  the  original  signals  is  to  remain  separable 
in   the   time  domain,  the  bandwidth  (BW)  of  the  system  must  be  large 
enough  to  allow  much  of  the  multiplexed  frequencies  to  pass.   As  an 
example  of  the  process  which  occurs  when  the  BW  is  limited,  let  us 
assume  that  the  audio  filter  allows  three  times  the  initial  BW  to  pass 
so  that  the  frequency  domain  result  is  as  shown  in  figure  6-4.     This 
can   be   inverted   to   produce   the   equivalent   time   domain   functions 
which  would   multiply   the  original  signals.     These  would  appear  as 
shown   in   figure   6-5. 
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Figure     6-4.         A     frequency     truncated     version     of     the     switched 
bandlimited    signal. 
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Figure   6-5.      Equivalent   time   domain   signals   which   are   obtained   by 
the   inverse   Fourier   transform   of  the   spectrum   in   figure   6-4. 


The  sample  at  tj  would   produce    1.14   Si(ti)  -  0.14  S2(M)  and  the 
sample  at  t2   would  produce   1.14  S2(t2)  -  0.14  Si(t2).     Clearly,  this  is 
a  significant  effect  and  the  exact  timing  of  the  sample  compared  to 
the   switching   function   is   also   quite   important. 

The  ideal  is  to  amplify  and  filter  at  RF,  then  multiplex  the 
signals,   finally   sampling   at   the   appropriate   rate   without   the   audio 
filter  in  place.     Note  that  the  number  of  points  per  image  decreases 
by  a  divisor  of  N  unless  the  number  of  points  is  changed  in  the  pulse 
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sequence.         The  sampling  rate  and  number  of  points  could  be 
controlled  by   selecting  a  specialized  pulse  sequence. 

An  alternative  method  of  obtaining  TDM   that  would  be 
effective  for  some  types  of  imaging  would  be  simply  to  delay  each 
successive  signal  until  the  prior  signal  had  decayed  to  noise.     The 
result  for  N  coils  would  be  an  acquisition  time  of  N  times  as  long. 
The  drawbacks  of  this  concept  include  the  difficulty  of  obtaining 
dispersionless  delays  of  several  microseconds  for  rf  signals   and  rapid 
pulse   sequences   in   which   acquisition   is   taking   place   nearly 
constantly. 

Frequency    Domain    Multiplexing 

The  diagram  of  TDM  with  limited  BW  brings  up  the  second 
implementation   discussed.      This   is   frequency   domain   multiplexing 
(FDM).     The  idea  is  to  simply  place  the  images  which  would  be 
produced  by   several   coils  into   adjacent  frequency  ranges.      This 
operation  can  be  done  most  cleanly  with  a  double  balanced  mixer 
and  a  low  frequency  sinusoid  as  the  local  oscillator.  As  an  example, 
the   simplest  implementation  for  two   signals   is   shown   in   figure  6-6. 
The   frequency   domain   picture   of  this   situation   looks   substantially 
like  the  filtered  TDM  except  now  the  separation  of  the  two  signals  is 
effected  by  their  frequency  ranges.  Figure  6-7  below  shows  the 

image  from  coil    1    unchanged   at  its  original  frequency  while  the 
signal  from  coil  2  is  replicated  at  higher  and  lower  frequencies.     If  a 
bandpass  filter  is  used  both  before   and  after  mixing   we  could  get 
three   images   in   this   same   frequency   range. 
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Figure  6-6.  A  simple  method  of  frequency  shifting  one  of  the 
bandlimited  signals  to  mix  the  signals  in  a  single  band  and  allow 
separation   by   frequency   position. 
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Figure   6-7.      Magnitude   Fourier   spectrum   of  the   two  frequency 
multiplexed    signals. 


An  alternative   method  of  performing  FDM  is  utilizing  an  extra 
set  of  gradient  type  windings   which  would  allow   for  a  stepped  static 
magnetic  field  positionally  linked  to  the  position  of  an  array  of  coils. 
The  end  result  is  that  the  signals  induced  in  each  coil  would  be 
mostly    separable    simply   by   the   frequencies   they   should   contain. 
Clearly  this  would  be  fairly  difficult  to  implement  since  it  would 
require   a   stairstep   magnetic   field.      Furthermore,   each   coil   would 
receive   signal   from   adjacent  regions   which   would   be   eliminated 
when   each   signal   was  filtered  to   its  expected  band. 
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Two   techniques   have  been   presented   for   simply   mixing   two 
signals  into  a  single  channel  while  allowing  separation  after 
acquisition.     The  fundamental  idea  is  to  orthogonalize  the  two  signals 
in  some  way,  sum  them  and  adjust  the  channel  to  permit  the 
increase  of  information.       From  chapter  4  it  is  known  that  errors 
resulting  in  coupling  of  the  information  can  be  corrected  if  the 
coupling    has    been    characterized. 

RF   Prefiltering 

There  is  also  a  simple  way  of  producing  the  same  kind  of 
effects  in  a  specialized  manner,  in  which  the  coils  are  matched  to 
narrowband   filters   and    are    arranged   in   the   frequency    encode 
direction.     The  effect  is  to  filter  noise  from  the  region  in  which  it  is 
unwanted  (where  there  is  no  signal).  Consider  figure  6-8  which 
shows   a  hypothetical   image  of  two   separate   samples.      The   frequency 
encode  gradient  is   as   shown,  resulting  in  the   separation  of  the  two 
signals   in  frequency.      An  example  line  in  the  frequency  direction 
through   the   image  is   shown  in   figure   6-9   with   a   representative 
separation  into  two  constituent  images  which  come  from  each  coil.     It 
is  clear  that  elimination  of  the  noise  in  the  band  of  the  other  coils 
signal,  as  shown  in  figure  6-10,  would  lead  to  an  increase  in  SNR  by 
V2   for  each  sample  region.     This  can  be  accomplished  with  filters 
which   are   offset   in   frequency   and   have   extremely   steep   skirts. 
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Figure   6-8.      A  hypothetical  image  of  two  separate   samples  from  two 
separate   coils,   with   the   frequency   direction   left   to   right 


To   demonstrate   the   efficacy  of  this   method,   images   were  made 
with  and  without  filters.     The  test  setup  consisted  of  four  loops  with 
decoupling   matching   circuits   designed   to   increase    isolation   between 
loops.      Four  low   impedance   preamplifiers   were   connected   followed 
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Figure  6-9.  Representative  line  in  frequency  direction  from  figure  6-8, 
with  separation  into  two  signals,  one  from  each  coil 
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Figure  6-10.  Representative  line  from  improved  image  showing  reduction 
in  noise  due  to  elimination  of  out  of  band  noise  by  filtering 
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by  two  90°   combiners.   This   formed  a  left  quadrature  breast  coil   and 
a  right  quadrature  breast  coil  the  outputs  of  which  were  fed  into  a  0° 
combiner  with  and  without  crystal  filters  matched  to  the  coil 
geometry  and  image  size.     A  schematic  of  this  is  shown  in  figure  6- 
11.     The  measured  SNR  in  the  region  of  the  breasts  increased  by 
about  40%  as  is  expected  by  theory,  and  images  of  a  phantom 
showing  this   improvement  are  shown  in  figure  6-12  with  the  filtered 
acquisition   at   the   bottom. 

In   figure   6-13   are   the   characteristics   measured   by   a   network 
analyzer  of  the   complementary   filters   used   to   obtain   the   images   of 
fig  6-12.     The  shoulders  of  the  filters  attenuate  the  signal  as  well  as 
the  noise,  however  only  at  severe  attenuations  is  the  SNR  adversely 
affected  and  it  is  assumed  that  this  signal  is  either  unwanted  and/or 
much  smaller  than  the  signal  from  the  adjacent  coil  and  filter.     This 
implementation   requires   that   the   filter   matches   the   coil   with   the 
image  bandwidth.     This  implies   that  only  a  particular  image  size  and 
bandwidth   is   allowable  for  a  certain  filter  set.     To  demonstrate   this 
consider  the   images   in   figure   6-14. 

The  top  image  is  of  a  rectangular  loader  with  long  axis  in  the 
frequency  encode  direction  with  a  20  cm     field  of  view,  while  the 
second  is  an  identical  image  but  with  a  36  cm  FOV.    This 
demonstrates  that  the  filter  is  intrinsically  matched  to   the   image 
bandwidth   not  the   spatial  characteristics  of  the  image.     This  is 
another  example  of  applying   knowledge   of  the   signals   to  a   specific 
implementation  to   obtain   greater  performance.      In   this   case   it  is 
necessary  to  know   the   signal  acquiring  characteristics  of  the  coil  in  a 


107 


it- 

V > 


Bandpass 
Crystal  Filters 


0°  0 

Combiner 


Output  to  MR  receiver 
Figure  6-11.  Pre-filtered  quadrature  breast  coils. 
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Figure   6-12.      Unfiltered  bilateral  image   at  top,   filtered   at  bottom 
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Figure   6-13.      Network   analyzer   transmission   plot   for   both   filters 
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Figure  6-1 4.  A  filtered  surface  coil  image  with  a  20  cm 
FOV  at  top  and  36  cm  FOV  at  bottom 
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spatial  sense,  and  it  is  assumed  that  the  distribution  of  sample  does 
not  greatly  change  the  overall  signal  profile. 

Another  difficulty  with  this  method  is  that,  unlike  the  previous 
methods   of  this   chapter,    coupling   between   channels   means 
information  is  irrevocably  lost.      Unlike   any   of  the  previous  methods 
in  which  reconstruction  is   done  after  acquisition,   the   addition  of 
signals  is   done  in  real  time  and  the  filters  simply  eliminate  any 
signals  from  outside  their  bands.     Therefore  coupling  produces  a 
direct  loss  in  SNR.     It  is  very  important  for  this  case  then  to  minimize 
coupling. 

Another  suggestion  has  been  made  in  the  literature  (42)  to  use 
this  general     technique  for  a  longitudinal  spine  array.     This  has  the 
potential  to  increase  SNR  but  not  to  optimize  it.     The  images  in  figure 
6-15   demonstrate   image   intensity   profiles   for   three   coils   arranged   in 
the  direction  of  the  main  magnetic  field.     The  difficulties  which  must 
be   overcome   in   this   arrangement   are   the   significant   signals   past  the 
physical   limits   of  the   coils   and  therefore   that   the   attenuation   profiles 
of  the  coils  must  match  the  relevant  signals  and  that  the  phase 
characteristics  of  the  filters  in  the  fall-off  region  are  controlled.     The 
phase  difference   between  filters   must  be   constant  through   the  region 
of   overlap. 

Conclusion 

The  acquisition  of  signals  from  multiple  coils  produces  the 
ability  to  increase  the  SNR  over  single  coil  systems.     The  most  flexible 
method  for  constructing  a  single  image  from  the  various  coils   is  by 
orthogonalizing  the  signals  in  some  way.     This  can  be  done  by  simply 
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Figure  6-1 5.  Images  from  each  of  three  surface  coils 

oriented  for  zero  mutual  inductance  in  the  z  direction 
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reproducing   the   receiver  and   never   letting   signals   from   different 
coils  be  mingled.     However,  any  method  which  allows  tagging  of 
signals  so  that  they  remain  separable  is  applicable  to  multi-channel 
and  single  channel  systems.     It  is  logical  to  expect  that  array  systems 
in   the  future   will   have   multiple   receivers   with   the   ability   to   increase 
the  effective  capacity  of  each  receiver.     This  may  be  done  by  any  of 
the   methods   described   above. 


CHAPTER  7 

MISCELLANEOUS  HARDWARE  IMPLEMENTATIONS 

Impedance    Matching 

Coils  for  NMR  and  MRI  use  have  historically  been  matched  to 
an  impedance  of  50  ohms.     This  is  one  of  the  standard  impedances 
for  test  equipment  and  amplifiers.     The  topics  in  this  dissertation  all 
deal   with   the   system   in   reception   and  therefore   impedance   matching 
as  power  matching   is  irrelevant.      It  is   easy   to   demonstrate  that 
maximum  power  is  delivered  to  a  load  for  which  the  impedance  is 
the  complex  conjugate  of  the  source  impedance.     Optimal   power 
matching  does   not  necessarily  mean  the  best  impedance  for 
reception. 

There  are  two  mechanisms  for  SNR  loss  produced  by  an 
impedance  mismatch.     The  first  is   that  the  preamplifier  has  a  noise 
figure  dependence  on  the  source  impedance.     The   second  is  an 
increase  in  transmission  line  loss  due  to  VSWR  effects.     The  greater 
the   impedance   mismatch    the    greater   the    standing    wave   voltages 
which  appear  across  some  small  resistances  in  the  cable. 

Many   matching   schemes   have  been   introduced  into  the 
literature   (20).      The   precise   method   is    substantially   unimportant,    but 
there   are   several   considerations   which   are   relevant.      First,   the 
matching  scheme  is  assumed  to  be  lossless.     This  is  assured  by 
allowing  only  reactive  elements   in   the   matching  network.      Second,  it 
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is   important   that   balanced   matching   techniques   are   used   (19).      These 
methods   serve   to  reduce  dielectric   losses   from  the   sample,   and   to 
reduce  the  susceptibility  to  external  noise  sources  (20). 

A   graphical  method  for  determining  the  component  values  for 
matching  has  been  derived  and   is  shown  below  as   an  example  of  the 
matching  procedure.     The   fundamental  part  of  the  probe   is   the   loop 
on  which  voltage  is  induced.     This  will  have  inductance  and 
resistance  and  is  therefore   shown   schematically  in  figure   7-1. 
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Figure  7-1.         A  circuit  model  of  the  fundamental  part  of  a  NMR  probe 

This   impedance   is   measured   at  the  frequency   at  which   the   match   is 
required   and  is   given   the   symbol   mZ0  .     The  procedure  for  the 
simple    matching    algorithm   using   two   reactive    components    configured 
schematically  as  shown  in  figure  7-2  is  given  below. 

1.  On  a  polar  graph,  as  shown  in  figure  7-3,  plot  the  vector 

rn'Z-0. 

2.  Draw  a  vertical  line  (lj)  through  the  head  of  this  vector. 


lm 
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Figure  7-2.     A  circuit  model  of  the  fundamental  part  of  a  NMR  probe 
with  reactances  necessary  for  matching  to  50  Q 
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3.  Draw  a  vertical  line  (12)  through  R0,  the  desired  end 

resistance  (usually   50  il). 

Icos  d) 

4.  Find  the  intersection  of  a  circle  of  radius      — —  (Ci)  with  l*. 

V  Rom 

The  distance  (shown  highlighted)  from  this  point  to  the 
end  of  the  vector   m"'Z-0  is  ^"'1. 

5.  Find  the  intersection  of  a  circle  of  radius   J—2—  (C2)  with  h. 

y  cosf 

The  distance  (shown  highlighted)  from  this  point  to  the 

real  axis  is  |;tm|. 
To   use   the  procedure   to   derive   values   by   measurement   instead   of 
calculation,  two  scales  should  be  used  on  the  graph-one  for  ohm's 
and  one  for  mho's.     The  two  sets  of  intersections  correspond  to  using 
a  capacitor  for   the   series   component   and   parallel   component,   while 
the  other  would  use   an   inductor  for  the   series  component  and  a 
different   capacitor   for   the   parallel   component. 
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Figure   7-3    A   graphical   procedure   for   deriving   reactances   necessary 
to  match  a  resistor  and  inductor  to  a  given  real  impedance  for  a 
specific   frequency    using   two   reactances. 


Roemer  et  al.  (39)  introduced  the  method  of  utilizing  the 
matching  network  as  a  method  of  producing  decoupling  between 
coils.     This  method  is  included  in  the  description  above  except  that 
the   original   source   impedance   has   had   some   inductive   reactance 
canceled   by   capacitive   reactance,   with   the   result   that  the   series 
inductance   has   the    same   magnitude   reactance   as    the   parallel 
capacitor.     This  works  well  for  elements  which  are  linear  coils.     The 
basic  concept,  as  previously  described,  is  to  make  the  effective  Q  of 
the  coil  low  by  loading  it  with  the  preamplifier  impedance.     This  in 
turn  reduces  the  ability  of  the  coil  to  serve  as  a  secondary  source  to 
another  receiving  coil.     This  method  is  most  effectively  used  with 
preamplifiers   in   close   proximity   to   the   coils;   otherwise   standing 
waves  can  couple  one  cable  to  another,  reducing  the  isolation  which 
renders    the    procedure    ineffective. 
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Impedance    preserving    combiners 

One  of  the   fundamental   tasks   discussed   in   the   previous   sections 
was  the  combination  of  signals  into  a  single  stream,   while 
maintaining  impedance  matching.     Although  this  is  basic  filter  design 
in  theory,  there  has  not  been  a  published  description  of  discrete 
component   phase    combiners    which    preserve    impedance    matching    in 
the   NMR   literature.      Several   different   implementations    are   discussed 
below. 

The    implementation    of   generalized    quadrature    requires 
combining  circuits  which  are  new  to  the  MRI  literature.     However, 
the   needs   are  very   similar  to  those   of     antenna   array  combinations 
for  ham  radio  operations.     The  general  philosophy   of  combiner  design 
is  that     one  must  produce  currents  which  differ  by  the  complex 
factor  desired.      For  equivalent  coils,   this   simply  requires  that  the 
coils   have   impedances   which   exactly   correspond   to   the   amplitudes 
and  phases  desired.     The  classic  case  was  described  by  Chen  et  al  in 
1983   (32),   and   produces   a   ninety   degree   phase   shift  between   two 
coils  which  are  placed  in  parallel  with  each  other.     One  coil  is 
matched  to  50  +J50  ohms  and  the  other  to  50  -  j50  ohms.     The 
parallel   equivalent  impedance   is   50   ohms   and  the  currents  produced 
are  45  -  -45  degrees  different  for  equivalent  coils.     This  can  be  easily 
generalized   to   any   phase   angle   and   splitting   amplitude,   but  it  may  be 
necessary  to  rematch  the  final  coil  to  approximately  50  ohms  if  the 
unit  coils  are  first  matched  to  50  ohms. 

The    straightforward    approach    mentioned    above   for   producing    a 
general   phase   combination   (with   equal    weighting)   is   described 
below,   except  the  original  coil  resistances  are  allowed  to  vary  from 
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50  ohms..     To  carry  out  this  method,  two  equations  must  be 
simultaneously   solved.     These   are    (z/'+Z.*"1)    =50  and  tan"1  — t   =  0  , 

where  $  is  the  angle  desired  for  combination.     These  equations  are 
solved  for  Z0  to  give  matching  for  each  coil  to  produce  the  desired 
effect.     The  drawback  of  this  method  is  that  there  is   a  requirement  to 
match  to  individual  coils  to  non-standard  impedances.     Nearly  all  test 
equipment  is  set  to  a  specific  nominal  frequency,  e.g.  50  ohms.     This 
makes    the    procedure    somewhat    tedious. 

In  addition,  it  is  possible  to  control  the  signal  weighting  in  a 
simple   manner  by   changing   the   reactances   of  the   combiner  circuit, 
by  increasing  one  and  decreasing  the  other,  it  is  possible  to  maintain 
the  phase  shift  of  interest.     This  results  in  fairly  small  changes  in  the 
final    impedance   but   can   produce    substantial    weighting    differences. 

Another  method  is  to  match  each  coil  independently  to  50 
ohms  and  then  use  a  impedance  preserving  phase   shifts  to  preset  the 
phase  to  the  desired  value.     This  can  be  followed  by  the  90° 
combiner   described   above,   which   results   in   50   ohm   output 
impedance.      The   shifts   mentioned  above   can  be   simply   realized  with 
the  following  algorithm.     Consider  the  simple  T  network  shown  in 
figure  7-4   .  The  phase  shift  of  this  network  when   loaded  with 
resistance  R  on  both  sides  is  -2  tan"'  —  j   different   than   if  there   were 

no  network  there  at  all.     The     impedance  is  preserved  at  value  R  if 
R2 +(a>L.f  =  2  — .        As  stated  above  this  is   straightforward  filter  theory 
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Figure  7-4.     A  "T"  Network  which  preserves  impedance  at  R  ohms, 
but  produces  a  phase  shift  compared  to  removal  of  the  network 

and  a  tc    network  could   be   easily   substituted. 


Transmission    synchronized    shielding 
Faraday  shields  have  commonly  been  used  in  MRI  in  recent 
years(17,  65,  66).     They  are  designed  to  allow     fairly  low  frequency 
gradient   magnetic    field    signals    to   penetrate    unattenuated   but   to 
prevent  rf  energy   from  passing  through   to   a  significant  extent.      They 
have  been   used  to  prevent  coupling  of  the  rf  coil  to  internal   gradient 
sets  and  to  some  parts  of  the  tissue  of  interest.     Recently  rf  shields 
have   been   used   as   rf  current  paths   to   decrease   radiation   of  energy 
by  effectively  making  a  dipole  into   a  quadrupole.     The   theory 
surrounding   the   use   of   shields   has   not   been   adequately   presented 
and    normally   simple   arguments   are   made   for   their   functions. 

The  rf  shield  should  be  a  physical,  frequency  selective  filter.     It 
must  be  a  low  pass  filter  to  allow  gradient  signals  to  penetrate  while 
attenuating   higher   frequency   magnetic    fields.      The   cutoff  frequency 
is  usually   somewhat  arbitrary   as   are   the   falloff  characteristics.      One 
common   method   of  producing  this  behavior  is   by   selecting   the 
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thickness  of  the  conductor  to  be  about  4  or  5  skin  depths  at  the  rf 
frequency   of  interest.      This   allows   a   nearly  unobstructed   image 
current  to  be  induced  at  rf  which  tends  to  cancel  the  field  on  the 
non-source  side  of  the  shield,  whereas  the  shielding  will  be 
ineffective  at  lower  frequencies.  Another  method  often 

incorporates   the   skin   depth   with   capacitive   bridges   between   a   ladder 
of  overlapping  plates  as  shown  in  figure  7-5.     This  provides   a  further 
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Figure  7-5.      Overlapping  fingers   of  copper  separated  by   a  dielectric 


dc  pole  for  the  description  of  the  effective  impedance  of  the  shield. 

A   new   shielding   technique   was   implemented   to   solve   a 
problem  of  rf  coupling  between  a  pair  of  quadrature  coils.     It  was 
desired  to  make   a  switchable  quadrature  coil   for  imaging   of  the 
human  breast.     This  was   a  challenge  because  each   quadrature  coil 
must  perform  in   two   different  environments.      Each   had   to   function 
when    the  other  is  either  on  or  off.     The  close  proximity  of  the  two 
coils  makes  this  difficult  because  of  the  coupling  of  energy  when 
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both  coils  are  on.     An  effective  rf  shield  was  placed  between  the  coils 
and   this   appears  to   solve  the  problem.      However,   there  was   another 
difficulty.       The  coils  were  receive-only  and  therefore  the  excitation 
of  the  spins  was  produced  by  the  surrounding  body  coil.     A  tip  angle 
artifact  appeared  in  the  image  in  the  vicinity  of  the  shield  because 
the  rf  of  the  body  coil  could  not  penetrate  the  shield.     This  problem 
was  eliminated  by  the  development  of  an  actively  controlled  shield. 
A  switching   voltage  is  available   which  is   synchronized  with   the 
transmitter   pulse   driving   the   body   coil.      A   capacitively   interdigitated 
shield     was  used  as  shown  above  in  figure     7-5.     It  was  observed  that 
interruption   of  the   strip   connecting   the   fingers   had   profound 
significance  for  the  rf  penetration  through  the  shield.     It  is  believed 
that   breaking    the    allowable   image   currents   into   areas    approximately 
half  the  size  of  the  source  loop  is  the  cause  for  this  loss  of  ability  to 
shield.     One  diode  circuit  was  inserted  at  a  break  in  the  connector 
strip  which  produced  an  effective  short  at  rf  when  the  diode  was 
reversed   biased   and   nonconducting   and   an  open  at  rf  when   the   diode 
was  forward  biased  and  therefore  conducting.     The  schematic  in 
figure     7-6  is  a  representation  of  this  and  figure  7-7     shows  the 
active  shield  in  operation  on  a   1.5T  Signa  scanner  with  the  left  image 
produced   with   the  shield   shorted  and   the   right  with   active   control. 
The  shield  was  in  contact  with  a  thin-walled  plastic  bottle  for  the 
images   of  figure  7-7. 
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Birdcage  Surface  Coil 
As   discussed  in  previous   sections,   considerable   interest 
remains  in  optimizing  surface  coils  for  MRI.     The  birdcage  is 
recognized  as  one  of  the  most  effective  volume  coils  ever  produced 


Figure  7-6.     A  schematic  representation  of  the  active  control  of  the 
gap    impedance 


and  although  attempts  have  been  made  to  produce  a  surface  coil 
which   mimics   birdcage   behavior   (67)   they   have   been   largely 
unsuccessful.     The  difficulty   seems  to  be   that  the  continuous  phase 
change   of  currents   in   a  transverse   direction   are   difficult  to   produce 
because   the   model   of  the   structure   as   a   transmission   line  requires   a 
loop  so  that  the  currents  circulate  in  the  structure.     A  segment  of 
transmission   line   has   a   fixed   standing   wave   if  terminated   on   either 
end. 
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Figure  7-7.  Surface  coil  images  of  a  bottle  -  actively 
controlled  shield  on  right,  no  control  on  left 
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A  solution  to  this  problem  has  been  found.     Since  the  birdcage  is  a 
discrete  model   of  a  transmission  line,   then  one  can  substitute   a  real 
transmission  line  for  part  of  the  structure  without  changing  its 
behavior  in  a  significant  way.     Consider  figure  7-8     which  depicts  a 
possible   implementation   of   this    concept. 


Figure  7-8.     A  Schematic  representation  of  a  birdcage  surface  coil 


The  structure  is  allowed  to  change  in  physical  size,  but  the 
characteristic   impedance  is   maintained.      This   impedance   is   of  course 
chosen  equal  to  that  of  the  real  transmission  line.     The  total  phase 
shift  across  the  structure  is  arbitrary  as  is  the  number  of  rungs.     Both 
of  these  should  be  based  on  the  region  to  be  imaged. 

The  result  of  this  implementation  is  a  closer  approximation  of  a 
limited   current   sheet   with   continuously    changing   phase,    than   any 
previous   surface  coil.     The  part  of  the   structure   which   is   transmission 
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line  will  be  essentially  lossless,  so  that  this  should  closely 
approximate   the  current  ideal   configuration  of  a  surface   coil. 


CHAPTER  8 
CONCLUSIONS 

As  knowledge  of  the  MRI  receiving  system  has  grown,  so  has 
the  ability  to  increase  performance  by  increasing  the  specificity  of 
our  system.     This  process  began  with  the  advent  of  surface  coils 
which   were   matched   to   the   region   of  interest,   proceeded   to 
quadrature   volume   coils,   then   quadrature   surface   coils   and   finally   to 
array    coils    and    quadrature    arrays. 

The  most  important  conclusions  to  this  work  are   summarized 
below.       First,  a  novel  approach  to  matching  coil  size  to  region  size 
was  investigated.      The   method  was   to   allow   continuous   adjustment  of 
the  coil  while  retaining  frequency  tuning  by  careful  element  design. 
Second,   the   nature   of  noise   correlation   and   the  relationship   with 
crosstalk  were  refined.      It  was  proven   that  noise   correlation   can  exist 
between  coils,  even  in  the  absence  of  any  crosstalk.     It  was  further 
proven   that  the  potentially  lossy  effects   of  crosstalk  can  be 
eliminated   by   including  knowledge   of  the   degree   and   nature   of  the 
coupling  into  the  combination  algorithm.     Third,   the  concept  of 
quadrature   coils   was   generalized   to   include   an   arbitrary   number   of 
coils,   with   arbitrary  phases   and  weights.   Though   similar  in 
theoretical   development   to   the   phased   array   work   (39),    no   other 
researcher   has   published    any    work   describing    single    channel 
combinations  of  signals.     Fourth,     the  general  orthogonalization  of 


127 


12! 


multiple   signals   into   a  single   receiver  was  conceptualized.      Several 
particular  methods   for   accomplishing   this   were   described.      Finally, 
several   new   instrumentation   concepts  were  described.      These   include 
a  simplified  matching  algorithm,   a  generalization  of  phase  combiners, 
a  transmission  synchronized  shield,  and  a  birdcage  surface  coil.     The 
tools  developed  as  well  as  the  new  methods  described  are  expected 
to  continue   to  have  relevance  in  future  research  in  the   area  of 
receiver   and   receiver  coil   design. 

The  use  of  arrays  of  coils  is  certain  to  increase,  since  this 
approaches  most  closely  the  optimum  performance.      The  number  of 
coils  used     in  array  systems  will  probably  increase.     The  currently 
conceivable   limit   is   the   number   required   to   perform   a   complete 
inversion  of  the  data,   so  that  the  number  of  coils  would  approach 
several  hundred.     The  use  of  multi-coil  systems  can  be  utilized  with 
adaptive   algorithms  which  enable   selection  of  coils  to  be  used  in 
acquisition.     A  possible  solution  would  involve  a  reduction  in  the 
number  of  gradient  steps  for  phase  encoding.     This   has  been   termed 
reduced   k-space   sampling.      This   provides   the   ability   to   acquire 
images  in  times  which  are  an  order  of  magnitude  or  more  faster.     The 
SNR  for  a  given  time,  however,  will  not  increase  nearly  as  strong  as 
this.     The  number  of  coils  can  allow  spatial  encoding  which  alleviates 
the  burden  of  phase  encoding,   however  the  averaging  effect  on   SNR 
from  this  encoding  will  be  eliminated  with  no  real  change  in  SNR  for 
a   given   acquisition   time. 

In   addition,   superconducting   receiver   coils   (69)   should   also 
increase   in   popularity.   These   have   the   ability   to   make   dramatic 
improvements   in   low   to   mid-field   imaging   where   the   coil   conductor 
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losses  dominate  for  non-superconducting  coils.      In  higher  field 
systems,  they  open  a  new  door  to  coil  design  through  which  coils  not 
heavily  loaded  by  the  sample  can  nonetheless  be  effective.     Target 
field  type   approaches   (70)   which  result  in  complicated   current 
designs   producing   near  optimal   fields   can   be   utilized   without  paying 
the  penalty   of  overwhelming   conductive   losses.      Methods   which 
require  precise  steep  boundaries  on  coil  sensitivities  can  be  realized. 

Three  and  four  Tesla  systems  have  been  routinely  utilized  for 
human  studies  for  some  time  (71)  and  will  probably  have  some  place 
in  the   future   clinical  environment.     Coil  designs   at  these  frequencies 
must  consider  radiation  effects  to  be  successful.     Designs  which  were 
extremely   ineffective   at   lower   field   strengths   may   function   well   at 
these  field  strengths.      Finally,   at  these  frequencies  the  eddy  currents 
induced   in   the  patient  may  compare  to  the   source  currents  in   the 
coil,   which  leads  to  the  interesting  difficulty  of  considering  the 
sample   to  be  part  of  the  receiving  system.     This  surely  makes   things 
more    complicated,    but   may    permit   the   opportunity    for   unrivaled 
local   images   within  the  body. 

It  is  the  hope  and  belief  of  the  author  that  the  research 
described  in  the  previous  chapters  is  applicable  to  most  of  the  areas 
expected  to  be  of  future  significance   in  NMR  receiver  hardware 
development. 
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